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NUMERICAL  ANALYSIS  OF  SEMICONDUCTOR  DEVICES 
USING  MICROCOMPUTERS 

1.  INTRODUCTION 

The  development  of  semiconductor  electronics  can  be 
measured  in  terms  of  improvements  in  device  switching 
speed,  cut-off  frequency,  geometric  dimensions  and 
complexity.  The  application  of  new  materials,  advances  in 
material  preparation  and  quality,  and  the  creation  of  new 
fabrication  technologies  have  all  been  necessary  to  meet 
the  ever  increasing  demands  of  the  industry.  Today,  for 
example.  Gallium  arsenide  (GaAs)  field  effect  transistors 
(FETs)  with  0.5  micron  gate  lengths  are  commonplace,  and 
current  research  efforts  are  directed  toward  devices  with 
even  smaller  features. 

Greater  device  and  material  complexity  has  brought 
with  it  the  need  for  improved  analysis  procedures.  In 
many  cases,  closed  form  solutions  to  the  often  nonlinear 
systems  of  equations  formulated  to  represent  these  devices 
are  not  tractable.  Designers  have,  quite  naturally, 


turned  to  the  use  of  computers  to  predict  device 
characteristics  and  to  gain  insight  into  device  physics. 
There  has  grown  up  an  effort  in  numerical  analysis  and 
device  modeling  that  parallels  developments  in  the  rest  of 
the  field. 

The  problems  that  this  segment  of  the  device  community 
face  include  limitations  placed  on  the  size  and  complexity 
of  devices  that  can  be  studied  by  computer  memory  size, 
machine  speed,  and  cost  of  computing.  Issues  such  as  the 
selection  of  boundary  conditions,  the  convergence  of 
solutions,  their  accuracy,  nonlinear  optimization 
strategies  and  the  like  must  be  jdalt  with  on  a  case  by 
case  basis.  The  cost  of  simulations  can  be  prohibitive  or 
can  limit  the  ranges  of  parameters  that  can  be  considered. 

The  purpose  of  the  work  reported  here  is  to  develop  a 
low  cost  modeling  technique.  The  work  is  directed  toward 
the  use  of  small  computers  which  normally  have  severe 
memory  restrictions  relative  to  the  capacity  of  larger 
machines.  A  further  goal  is  to  employ  an  approach  where 
the  convergence  of  the  solutions  and  their  accuracy  can  be 
insured.  The  intention  is  to  provide  a  user  oriented 
package  that  can  be  readily  reconfigured  to  allow  for  the 
numerical  analysis  of  arbitrary  devices.  The  results  of 
these  efforts  are  reported  here. 

One  of  the  earliest  use.  *  '  mputers  to  solve  a 
semiconductor  device  problem  was  the  self-consistent 


iterative  scheme  for  one-dimensional  steady-state 
calculation  reported  by  Gummel  (1)*.  Two  years  later  a 
transient  numerical  analysis  for  solving  Gunn  diode 
problems  was  published  by  NcCumber  and  chynoweth  (2). 
From  1968  to  1973,  many  papers  reported  techniques  for  the 
one-dimensional  analysis  of  a  variety  of  devices  (3-7). 
Among  them  the  work  of  Hachtel, et.  al  (6)  gives,  perhaps, 
the  most  general  and  efficient  analysis  for  junction 
device  modeling. 

For  devices  with  comparatively  simple  geometries,  such 

as  bipolar  transistors,  one-dimensional  numerical  analysis 

might  be  adequate.  However,  modem  designs  involve 

devices  so  complicated  that  this  analysis  fails  to  offer  a 

sufficient  and  accurate  evaluation  of  device 

characteristics.  Therefore,  two-dimensional  or  even 

three-dimensional  numerical  analysis  becomes  necessary  for 

device  designers.  The  first  two-dimensional  solutions 

appeared  in  1969,  reported  by  Slotboom  (8)  and  Kennedey 

and  O'Brien  (9),  separately.  Following  these,  the  FET 

became  a  popular  device  for  two-dimensional  investigation 

(10-19).  In  1976,  Yamaguchi  and  Kodera  (20)  reported  a 

sucessful  treatment  to  reduce  the  computing  effort  by 

introducing  a  regional  model. 

*  The  number  in  parentheses  in  the  text  indicate 
references  in  the  bibliography. 


Up  to  this  point  the  numerical  techniques  used  in 
device  simulations  were  based  on  the  finite  difference 
method.  The  physical  model  used  almost  universally  is  the 
well  known  diffusion  model  which  includes  Poisson's 
equation  and  the  current  continuity  equation.  In  other 
words,  the  nature  of  the  problem  is  to  apply  a  finite 
difference  scheme  to  solve  the  system  of  nonlinear  partial 
differential  equations  as  a  boundary  value  problem. 
Seqential  iteration  is  applied  until  the  solution 
converges . 

The  finite  element  method  was  introduced  in  device 
modeling  literature  in  1974  (21-27)  and  now  provides  an 
alternative  to  the  finite  difference  method  for 
discretization.  Reasons  for  using  the  finite  element 
method  include  (1)  ease  in  treating  irregular  geometries, 
(2)  automatic  conservation  of  current,  (3)  ease  in  local 
refinement  of  mesh  size,  and  (4)  ease  in  constructing 
higher  order  approximations.  However,  this  is  not  to  say 
that  the  finite  difference  method  has  lost  its  value  in 
device  modeling.  Rather,  with  its  simplicity  and 
reliability,  this  method  is  still  useful  in  most  cases. 

A  common  assumption  made  by  many  authors  is  that  the 
device  under  consideration  can  be  described  by  a 
two-dimensional  model.  This  is  generally  true  and  the 
result  is  adequate  for  engineering  purposes.  However  for 
those  devices  where  the  cross-section  cannot  be  taken 
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uniform  over  the  entire  third  dimension,  three-dimensional 
analysis  becomes  a  necessity.  Up  to  the  present  time, 
this  kind  of  treatment  is  rare  in  the  literature  (28). 
Recently,  a  program  named  FIELDAY,  developed  by  IBM,  was 
reported  (26).  This  program  simulates  semiconductor 
devices  of  arbitrary  shape  in  one,  two  or  three  dimensions 
operating  under  transient  or  steady-state  conditions. 
Also,  this  program  enables  users  to  rapidly  generate  new 
models  and  analyze  the  results.  It  is  believed  to  be  the 
most  general  and  efficient  algorithm  for  device  design 
purposes  at  the  present  time. 

Computer  storage  requirements  and  processing  costs 
have  restricted  computer-aided  design  activity  to  those 
institutions  capable  of  bearing  the  expense.  However  the 
recent  introduction  of  personal  computers  has  made  it 
feasible  to  do  certain  classes  of  device  simulations  very 
cost  effectively.  To  achieve  this  purpose,  a  simulation 
technique  must  be  found  that  fits  with  the  comparatively 
small  CPU  memory  space  available  for  calculation.  A  newly 
developed  modified  decoupled  method  has  been  successful  in 
satisfying  this  requirement.  Based  on  this,  a  user 
oriented  program  suitable  for  the  simulation  of 
metal-semiconductor  devices  has  been  developed  and  is 
reported  here.  The  approach  can  be  applied  to  any 
semiconductor  device. 


In  this  algorithm,  a  five  point  finite  difference 
scheme  was  chosen  for  simplicity.  It  can  be  shown  that 
the  finite  element  method  can  also  be  applied  equally 
well.  The  Newton-Raphson  iteration  scheme  combined  with  a 
direct  matrix  solution  technique  will  always  yield  a 
converged  result,  especially  for  diagonally  dominant 
matrix  systems  of  the  type  encountered  here.  Of  course, 
there  are  some  natural  limitations  which  dominate  the 
accuracy  of  the  result  and  the  device  area  that  can  be 
simulated,  such  as  the  round-off  error,  the  mesh  Debye 
length  limitation,  the  size  of  the  memory  available,  the 
accuracy  of  boundary  values,  and  so  on. 

An  additional  purpose  of  this  work  is  to  develop  a 
new  family  of  devices,  namely  the  geometrically  controlled 
metal-semiconductor  diode.  With  the  aid  of  the  program 
mentioned  above,  two  new  types  of  diodes  have  been 
investigated. 

1.1  NUMERICAL  ANALYSIS  (29) 

The  diffusion  model  is  generally  accepted  in  device 
modeling  and  requires  no  further  discussion.  To  obtain  a 
solution,  Poission's  equation,  the  current  equation,  and 
the  continuity  equation  are  solved  together  numerically, 
subject  to  the  appropriate  boundary  conditions.  In  this 
work,  for  simplicity,  generation  and  recombination 
processes  are  neglected.  Likewise,  the  formulation  is 
restricted  to  a  static  analysis  and  so  the  time  dependent 


terms  are  also  omitted.  Therefore,  the  continuity 
equation  can  be  combined  with  the  current  equation, 
thereby  reducing  the  problem  statement  to  two  equations  in 
two  unknowns;  the  potential  and  carrier  concentration  at 
each  point  in  the  device. 

There  are  three  steps  used  in  numerical  analysis:  (1) 
discretization,  (2)  solution  and  (3)  iteration. 
Discretization  describes  the  process  of  converting  a 
continuum  problem  into  its  discrete  version  (30).  For 
this  purpose,  either  the  finite  difference  or  finite 
element  method  can  be  chosen.  In  either  case,  two  sets  of 
nonlinear  algebraic  equations  result  after  discretization, 
one  from  Poisson's  equation,  the  other  one  from  the 
current  equation.  Step  two  is  to  solve  the  above  equation 
sets.  Two  different  methods  are  accepted,  the  coupled 
method  and  the  decoupled  method.  With  the  coupled  method 
the  two  sets  of  equations  are  solved  simultaneously.  In 
the  decoupled  method  they  are  solved  serially.  For 
iteration,  either  the  direct  iteration  or  the 
Newton-Raphson  method  is  available.  This  will  be  discussed 
further  in  chapter  3. 

The  advantages  and  disadvantages  of  the  finite 
difference  and  finite  element  methods  have  been  metioned 
previously.  In  this  work,  the  finite  difference  method  has 
been  chosen  to  simplify  the  programming  work  and  to  allow 
more  attention  to  be  paid  to  developing  a  new  approach  to 


be  described  later.  The  coupled  method  requires 
approximately  4  to  9  times  the  amount  of  memory  needed  by 
the  decoupled  method,  for  unipolar  and  bipolar  devices, 
respectively,  and  is  good  for  strongly  nonlinear  types  of 
problems.  The  decoupled  method,  which  solves  weakly 
nonlinear  problems  sucessfully,  is  the  most  widely  used 
approach.  Direct  iteration  is  simple  but  has  only  a 
first-order  rate  of  convergence,  while  the  Newton-Raphson 
method  requires  a  more  complicated  procedure  but  premises 
second-order  convergence.  Since  numerical  analysis  using  a 
minicomputer  is  emphasised  in  this  work,  special 
accomodation  is  needed  for  the  limited  memory  space 
available.  Considering  the  massive  memory  required  by  the 
coupled  method,  we  have  rejected  it  in  favor  of  the 
decoupled  method.  However,  with  this  method,  convergence 
is  not  guaranteed. 

In  chapter  3,  a  new  treatment  to  manipulate  and  solve 
this  problem  is  reported,  called  the  modified  decoupled 
method.  This  method  employs  the  Newton-Raphson  technique 
and  so  it  is  suitable  for  solving  weakly  or  strongly 
nonlinear  problems.  The  results  are  as  good  as  those 
obtained  using  the  coupled  method. 

In  the  numerical  analysis  of  semiconductor  devices 
there  are  some  additional  matters  to  be  considered.  Since 
device  simulation  is,  basically,  the  solution  to  a 
boundary  value  problem,  the  accuracy  of  the  boundary 


values  used,  such  as  carrier  concentration  and  potential 
at  the  contacts,  directly  affects  the  accuracy  of  the 
results.  Unfortunately,  precise  boundary  conditions  are 
not  well  understood  yet.  For  example,  the  exact  carrier 
concentration  at  a  Schottky-barrier  contact  is  not  known. 
Therefore,  contact  problems  need  to  be  studied  in  more 
detail.  The  same  argument  can  be  applied  to  the  assumed 
carrier  velocity-electric  field  and  diffusion 
coefficient-electric  field  relations  which  depend  on  the 
semiconductors  used. 

Mesh  spacing  is  another  factor  which  affects  the 
accuracy  of  the  numerical  results.  It  will  be  shown  in 
chapter  3  that,  in  general,  the  Debye  length  restriction 
is  a  good  reference  to  guarantee  convergence. 

Once  the  problem  has  been  formulated,  reasonable 
boundary  conditions  defined,  computer  programming 
completed  and  numerical  solutions  obtained,  one  is  faced 
with  the  task  of  verifying  the  validity  of  the  results. 
There  are  three  means  available  to  the  device  designer  to 
check  the  self-consistency  of  the  results  of  the 
simulation. 


1.  Conservation  of  current:  The  current  into  the  device 
should  be  equal  to  the  current  out  of  the  device. 

2.  Comparison  with  the  one-dimensional  result: 
Sometimes  one-dimensional  results  are  well  known 
analytically  and/or  phsically.  These  are  good 
references  for  comparison. 

3.  Symmetry:  With  symmetrical  or  anti  symmetrical 

geometries  the  results  should  conserve  the 
corresponding  symmetry  properties  which  can  be 
used  as  a  reference  as  to  the  self-consistency  of 
the  solution. 

1.2  APPLICATIONS  IN  SEMI  CONDUCTOR  DEVICES 

In  this  work,  a  two-dimensional  self-consistent 
computer  program  has  been  written  for  a  64K-byte  HP  9825B 
desktop  calculator  to  investigate  the  electrical 
properties  of  devices  with  arbitrary  contact 
configurations  and  general  boundary  conditions.  It  has 
been  applied  to  two  new  types  of  devices  which  show 
promise  as  microwave  and/or  millimeter-wave,  (mm-wave), 
devices . 

GaAs,  Schottky-barrier  diodes  (SBDs)  are  widely  used 
in  microwave  and  mm-wave  applications  such  as  detectors, 
mixers,  limiters  and  multipliers.  Especially  for  mm-wave 
applications,  usable  devices  require  a  low  noise  figure, 
low  parasitic  spreading  resistance  and  low  barrier 
capacitance.  GaAs  is  a  material  which  lends  itself  to 


meeting  ell  these  requirements  (31).  The  one-dimensional 
behavior  of  metal-semiconductor  rectifier  systems,  i.e. 
Schottky-barriers  are  well  understood.  However  in  many 
devices,  such  as  the  metal-semiconductor  field  effect 
transistor  (MESFET),  the  contact  geometries  are 
distributed  and  so  a  two-dimensional  model  for  the  device 
is  needed.  Moreover,  a  model  for  GaAs  devices  must 
accomodate  nonlinear  material  properties,  such  as 
field-dependent  velocity  and  the  diffusion.  Although  it 
may  be  possible  to  formulate  the  general  problem 
analytically,  the  solution  of  the  resultant  nonlinear  set 
of  equations  can  best  be  obtained  numerically. 

Mixer  diodes  employing  Schottky-barrier  contacts  are 
found  to  be  particularly  useful  in  high  resolution  radar, 
microwave  communications,  and  radio  astronomy.  Typically 
their  cutoff  frequencies  are  in  the  range  of  400  to  1000 
GHz.  Some  of  these  diodes,  such  as  beam  leaded  ones,  can 
go  beyond  1000  GHz  (32).  In  most  of  the  cases,  low 
turn-on  voltage  is  desired  to  increase  detector 
sensitivities  in  some  laboratory  test  equipment  and 
especially  in  systems  where  local  oscillator  power  is 
limited  (33).  Here,  the  turn-on  voltage  is  defined  as  the 
applied  voltage  needed  to  obtain  a  particular  forward 
current,  for  example  1  yA  (Figure  1.1).  The  turn-on 
voltage  is  closely  related  to  the  barrier  height  of  the 
rectifying  system.  For  low  turn-on  voltage,  low  barrier 
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height  is  required.  Uhder  regular  fabrication  procedures, 
for  a  given  Schottky  metal  and  a  particular  semiconductor 
material,  the  barrier  height  should  be  a  fixed  value. 
However  the  number  of  choices  for  the  metal  is  limited  by 
technical  considerations.  Some  barrier  height  lowering 
techniques  can  be  applied  to  meet  this  requirement,  such 
as  high  temperature  processing  (34),  interlayer  charge 
control  using  ion  implantation  (35),  and  the  planar  doped 
barrier  (36). 

The  work  presented  here  reports  a  new  method  to 
control  the  tum-on  voltage  of  rectifying  contacts  using 
conventional  Schottky  metal  barriers  on  arbitrary 
semiconductors.  This  new  approach  involves  including  the 
contact  geometry  as  part  of  the  design.  This  method  has 
been  proven  to  be  simple,  effective  and  can  be  achieved  at 
low  cost. 

1.2.1.  Gap  Controlled  Metal-Semiconductor  Diode 

The  first  device  under  consideration  is  called  the  Gap 
diode.  Its  current-voltage  characteristic  is  controlled  by 
a  geometric  gap  between  adjacent  Schottky  barrier  regions 
at  the  blocking  contact  on  one  side  of  the  diode. 
Physically,  it  can  be  expected  that  the  effective  barrier 
height  of  the  Gap  diode  can  be  adjusted  to  any  desired 
value  between  zero  and  the  characteristic  barrier  height 
of  the  Schottky  metal  on  the  semiconductor  material.  This 
means  that  the  applied  voltage  required  to  turn  on  the 


■k. 


diode  is  adjustable  simply  by  changing  the  width  of  the 
gap.  This  device  has  been  demonstrated  experimentally. 
Striking  similarity  between  the  predicted  and  experimental 
I-V  curves  was  observed  (37). 

1.2.2.  Antisymmetric  Metal-Semiconductor  Diode 

The  second  new  device  under  investigation  is  called 
the  Antisymmetic  diode.  Just  as  the  Gap  diode  is  a 
modified  version  of  the  Schottky-barrier  diode,  the 
Antisymmetric  diode  is  a  modified  version  of  an 
antiparallel  SBD  pair.  It  has  been  shown  (38)  that  the 
antiparallel  SBD  pair  can  be  pumped  subharmonically  at  one 
half  of  the  local  oscillator  frequency  and  can  be  used  as 
a  subharmonic  mixer.  In  microwave  and  mm-wave  frequency 
applications,  this  has  several  advantages  including 
relaxed  requirements  on  the  local  oscillator  frequency  and 
power,  a  reduced  noise  figure  and  an  improved  burn-out 
rating.  However  two  disadvantages  of  this  structure  are 
the  difficulty  in  fabricating  a  completely  balanced  diode 
pair  including  the  necessary  circuit  and  contact,  and 
secondly  the  turn-on  voltage  is  not  adjustable.  The 
antisymmetric  diode  described  here  is  one  approach  to 
solve  these  problems. 

1.3  SCOPS  OF  THIS  STUDY 

The  main  purposes  of  this  study  are  to  develop  a 
computer-aided  analysis  for  two-dimensional  metal- 
semiconductor  devices,  and  to  use  it  to  design  new  types 


of  microwave  and  mm-wave  devices.  In  Chapter  2 
one-dimensional  metal-semicondutor  physics  is  reviewed 
with  emphasis  on  the  depletion  approximation  and 
deviations  from  the  ideal  case.  Chapter  3  describes 
computer-aided  simulation  techniques  for  one  and  two 
dimensions.  A  modified  decoupled  method  is  developed  for 
the  solution  of  coupled  nonlinear  equations  suitable  for 
use  on  small  computers .  Some  technical  problems  are 
discussed,  such  as  boundary  values,  mesh  sizes,  and  so  on. 
A  one-dimensional  example  of  a  Schottky-barrier  junction 
is  presented  and  the  results  compared  with  those  of  the 
one-dimensional  depletion  approximation. 

The  analysis  method  is  applied  to  a  new  microwave 
diode  called  the  Gap  diode.  Chapter  4  presents  the 
results  of  the  simulation  and  compares  them  to 
experimental  results .  Chapter  5  reports  another  new 
diode,  the  antisymmetric  diode  which  is  another  potential 
microwave  and  mm-wave  device.  Theoretical  results  are 
presented. 


2.1  INTRODUCTION  (39) 

In  this  chapter,  the  characteristics  of  the  one- 
dimensional  metal-semiconductor  diode  are  discussed  in 
order  to  develop  a  basic  understanding  about  the  diode  and 
to  provide  background  for  two-dimensional  analysis. 


Metal-semiconductor  contacts 


are  also 


called 


Schottky-barrier  contacts.  In  the  area  under  the  Schottky 
metal  contact,  there  exists  a  region  almost  depleted  of 
free  carriers  compared  to  the  concentration  of  electrons 
supplied  by  donor  atoms.  An  electric  field  exists  in  this 
depletion  region  supported  by  the  fixed  charge  of  the 
ionized  immobile  donors.  The  extent  of  the  depleted 
region,  W,  under  zero  bias  is  called  the  depletion  length. 

Figure  2.1  (a)  shows  an  n-type  semiconductor  diode 
comprising  a  Schottky  metal  contact  and  an  ohmic  contact 
applied  to  the  semiconductor.  The  band  diagram  and  the 
electron  current  component  corresponding  to  forward,  zero 
and  reverse  bias  cases  are  shown  in  Figure  2.1  (b),  (c) 
and  (d),  respectively.  Here  the  diode  length,  L,  (the 
length  of  the  bulk  region),  is  assumed  to  be  greater  than 
the  depletion  length,  WQ.  An  ideal  ohmic  contact  is 
defined  as  a  thermal  equilibrium  contact  without  a 
potential  barrier  and  with  an  infinite  surface 
recombination  velocity  (40).  The  potential  difference 


between  the  Schottky  metal  contact  and  the  edge  of  the 
depletion  region,  denoted  by  V  ,  is  called  the  contact 
potential,  built-in  potential,  or  diffusion  potential.  In 
Figure  2.1,  the  barrier  height  d>w  is  the  difference 

O 

between  the  Schottky  metal  work  function  and  the  electron 
affinity  of  the  semiconductor.  Vf  and  Vr  are  the 
magnitudes  of  the  forward  and  reverse  bias  voltages, 
respectively.  The  elecron  current  from  the  semiconductor 
to  the  Schottky  metal  contact  is  denoted  by  and  that 
from  Schottky  metal  contact  to  semiconductor  is  denoted  by 
Jms  .  The  barrier  height  as  seen  from  the  metal  side 
stays  the  same  over  a  wide  range  of  applied  voltages,  V, 
between  VQ  and  reverse  breakdown  voltage.  So  Jms  remains 
constant.  On  the  other  side,  Jsm  which  is  the  current 
density  from  the  semiconductor  to  the  metal  must  equal 
Jms  when  no  bias  is  applied.  This  is  shown  in  Figure  2.1 
(c).  The  forward  bias  voltage  Vf  decreases  the 
diffusion  potential  and  allows  electrons  to  surmount  the 
barrier  from  the  semiconductor  to  the  metal,  so  that 
Jsm  increases.  The  reverse  bias  voltage  Vr  increases 
the  diffusion  potential  and  decreases  the  number  of 
electrons  able  to  surmount  the  barrier  from  the 


semiconductor  to  the  metal,  so  that  Jsm 

large  V„  ,  J  will  come  close  to  zero, 
r  sm 


decreases .  For 
In  this 


case,  J _  can  be  understood  to  be  the  saturation  current 

ms 

density. 


There  are  two  mechanisms  for  charge  transport  over  the 
potential  barrier  of  Figure  2.1.  If  electrons  have  enough 
thermal  energy  to  surmount  the  barrier  then  according  to 
thermionic-emission  theory  (39),  the  total  current 
density,  Jn  ,  is  related  to  the  applied  voltage,  V,  by 

Jn  =  A*T2exp(-qV0/kT) [exp (qV/kT) -1]  (2.1) 

*  *  _2  _  2  * 

where  A  =  120  (m  /m)  Amp  cm  °K  ,  for  which  m  /m  is 

the  ratio  of  semiconductor  effective  mass  to  free  electron 

mass.  For  GaAs,  this  ratio  is  0.072  (39).  The 

temperature  is  in  degrees  Kelvin.  According  to  diffusion 

theory  (39,41),  the  total  current  can  be  written  as 

Jn  *  qNdynEmaxexpt”qVc/kT)  (exP(<3v/kT) -1]  (2.2) 

where  g  is  electronic  charge,  Nd  is  the  doping 
concentration  of  the  bulk  region,  vn  is  the  electron 
mobility,  and  Emax  is  the  field  strength  at  the  Schottky 
contact.  These  two  completely  different  mechanisms 
predict  similar  results  and  both  are  accepted.  In  fact, 
it  is  possible  to  make  SBDs  in  which  the  I-V  curves  are 
very  close  to  these  predictions,  especially  in  the  forward 
bias  case  (39).  However,  effects  such  as 
quantum-mechanical  tunnelling  through  the  barrier  and  hot 
electrons,  to  name  a  few,  cause  departures  from  the 
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aforementioned  ideal  results.  This  will  be  explained 
later. 

In  the  practical  case,  under  forward  bias,  the  I-V 
curve  increases  exponentially  as  theory  predicts. 
Usually,  significant  current  conduction  will  be  detected 
beyond  the  tum-on  voltage  which  is  very  close  to  the 
built-in  voltage.  On  the  other  hand,  under  reverse  bias, 
the  I-V  curve  of  the  real  case  is  not  so  close  to  the 
predicted  results.  In  fact  the  reverse  current  increases 
with  the  voltage.  Eventually  the  avalanche  breakdown 
process  will  dominate  and  this  results  in  a  sharply 
increasing  current  when  V  approaches  -VBD  .  This  will  be 
dissussed  later.  A  typical  I-V  curve  for  SBD  is  shown  in 
Figure  2.2. 

2.2  ONE-DIMENSIONAL  APPROXIMATION 

The  usual  treatment  of  the  metal-semiconductor 
junction  considers  a  one-dimensional  configuration  in 
which  the  contact  metal  is  infinite  in  extent.  The 
potential  distribution  in  the  bulk  semiconductor  region 
can  be  calculated  from  Poisson's  equation  using  a  suitable 
approximation  called  the  depletion  approximation. 

The  potential  and  the  field  distribution  are  described 
in  one-dimension  by  Poisson's  equation: 

Hr  *  I  <n"Nd) 


(2.3) 


A  voltage  V  is  applied  to  the  diode  shown  in  Figure  2.3. 
If  at  x  =  0,  the  potential  is  zero,  then  the  potential  at 
W  will  be  Vo  -  V,  due  to  the  definition  of  W.  Beyond  the 
depletion  edge  x  *  W,  i.e.  outside  the  depletion  region, 
the  field  is  very  close  to  zero,  since  here  the  electron 
density,  n,  is  so  close  to  Nd  that  the  right  hand  side  of 
equation  2.3  is  effectively  zero,  and  the  potential  is 
constant.  Inside  the  depletion  region,  n  is  so  small  that 
its  value  can  be  neglected,  and  the  field  intensity 
increases  linearly  from  zero  at  the  edge  of  depletion 
region  to  a  maximum  value  at  the  Schottky  metal  contact. 
The  potential  here  can  be  described  by  a  quadratic 
function  of  x.  This  is  a  good  approximation,  since  for 
the  electrons  to  climb  the  potential  barrier,  q  (Vo  -  V), 
the  value  of  (n  -  Nd)  must  change  drastically  in  a  very 
short  distance  around  the  depletion  edge. 

Explicitly,  the  boundary  conditions  are 
Potential : 


v  -v 


W  <  x  <  L 


2 


carrier  concentration  : 


n  -  0 


0  <  X  <  W 


n  «  N . 


W  <  x  <  L 


Inside  and  outside  the  depleted  region  equation  2.3 
becomes : 


-!“d 


0  <  x  <  w 


(2.5 


W  <  X  <  L 


Using  the  following  relation  as  a  boundary  condition 


SI-- 


W  <  x  <  L 


(2.6 


yields 


S  -  -  f  "d'*-"’ 


o  <  x  <  w 


and  so 


<Mx) 


?  Nj (4  x2 


0  <  x  <  w 


In  order  to  satisfy  the  boundary  condition  at  x  *  W,  the 
value  of  V  should  be 


w 


2e (VQ-V) 
<3  Nd 


1/2 


(2.9) 


Equation  2.7  shows  that  the  field  strength,  E  *  -  df  /  dx, 
at  x  *  0,  i.e.  the  surface  of  semiconductor,  is  maximum. 


max 


I  "a" 


2qHd(V0-V) 


1/2 


2(Vq-V) 


(2.10) 


W 


and  decreases  linearly  in  x  to  a  minimum  value  zero  at  x  * 
W.  The  direction  of  the  field  is  from  the  semiconductor 
to  the  metal.  The  total  charge  in  the  depletion  region  is 


Q 


e  EmaxA 


qNdWA 


and  the  junction  capacitance  is 


(2.111 


c  - 


(2.12; 


where  A  is  the  diode  area. 

2.3  CURRENT-VOLTAGE  MODEL 

As  mentioned  previously,  there  are  two  approaches  to 
describe  the  I-V  relation  of  a  Schottky  diode,  i.e.  the 
thermionic-emission  and  the  diffusion  model.  Since  the 


diffusion  model  can  be  generalized  and  applied  to  a 
two-dinensional  structure,  a  detailed  study  of  it  is 
worthwhile.  The  thermionic-emission  model  is  omitted 
here. 

The  expression  for  current  in  terms  of  drift  and 
diffusion  components  is 


Jn  "  +  qDn  35 1 


(2.13 


where  positive  J  means  in  the  positive  x  direction,  D. 
is  the  diffusion  coefficient  for  electrons  and  is  taken 
independent  of  position.  The  problem  is  to  find  an 
expression  for  Jn  in  terms  of  the  applied  voltage  V. 
Before  starting  to  calculate  the  current,  it  should  be 
noted  that  the  diffusion  model  is  based  on  the  following 
assumptions : 

1.  The  current  through  the  barrier  is  due  to  the 
combined  effects  of  drift  and  diffusion  in  the 
depletion  region. 

2.  The  barrier  height  is  much  greater  than  kT/g. 

3.  The  carrier  concentrations  at  x  »  0  and  W  are 

unaffected  by  the  forward  bias  voltage. 

The  simplest  expression  relating  Dn  and  un  is  the 
Einstein  relationship. 
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Since  the  static  case  ia  assumed,  dJn  /dx  =  0,  and  Jn  is 
a  constant.  Rewriting  aquation  2.13  by  adding  a  factor 
exp(-  qi|>  At),  together  with  equation  2.14,  gives 

Jn  exp(-q^AT)  -  Dn  [n  exp(-qiJ>AT)  (2.15) 

When  V  =  0,  there  is  no  current  (  J  *  0  ).  At  x  *  W,  the 

n 

potential  is  the  built-in  value,  *  Vo  and  n  =  Nd  ;  at 
x  =  0,  ip  =  0,  and  n  =  n(0).  So  from  equation  2.15  one 

obtains, 

n(0)  -  Nd  exp  (-qvoAT)  (2.16) 

From  assumption  3,  equation  2.16  will  be  applied  in  all 
cases.  With  the  aid  of  equations  2.8,  2.9,  2.10,  and 

integrating  with  respect  to  x  from  0  to  W,  yields  equation 

2.2. 


2.4  EXACT  ANALYSIS 

Instead  of  the  depletion  approximation,  sometimes  an 
exact  analysis  is  desired,  especially  for  the 
two-dimensional  case.  The  charge  density  and  potential 
are  related  by  Poisson's  equation.  See  equation  2.3.  The 
current  density,  in  terms  of  drift  and  diffusion,  is 
described  in  equation  2.13.  If  the  Einstein  relationship 
is  assumed,  a  simple  expression  for  the  field  at  the 
contact,  i.e.  Emax  ,  can  be  found. 
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If  Jn  *  0,  such  as  is  obtained  for  a  bias  voltage 
less  than  the  turn-on  voltage,  or  in  reverse  bias, 
equation  2.13  can  be  rewritten  as 

gj  In  exp(-q^/kT) ]  *  0  (2.17) 

Integrating  and  applying  the  boundary  condition  at  x  =  L, 
n  *  Na  ,  where  the  potential  is  Vo  -  V,  one  obtains 

N(x)  *  Nd  exp[-  <-<|H-vo-V)]  (2.18) 

Substituting  n(x)  in  Poisson's  equation,  one  obtains, 
after  sone  Manipulation, 

3x  <S>2  ■  HT*  |*“ptir  «*-Vv>l-i|  g  (2-is) 


Integrate  this  equation  from  x  *  0  to  L  to  get 


(£i) 

'dx 


2 

x«L 


at 


.  kr 
e  q 


v  +  VQ) 


(2.20) 


For  a  Schottky  diode  in  which  the  length  of  the 

semiconductor  is  ouch  greater  than  the  depletion  length, 

i.e.  L  >>  W  . 

o 


A  I 


I 


2.2 
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the  value  of  E^*  will  be 


Fr4  <vv-  f  >] 


1/2 


(2.22) 


Comparing  this  with  the  expression  for  Emax  derived  from 

the  depletion  approximation,  equation  2.10,  the  only 

difference  is  the  term  kT/q.  Bowever,  from  assumption  2, 

kT/q  <  <  Vo,  so  E _  in  both  cases  is  identical,  when 

max 

equation  2.21  holds. 

Even  in  the  one-dimensional  case,  closed  form 
solutions  for  n  and  are  not  attainable,  because  of  the 
nonlinear  nature  of  equation  2.19.  On  the  other  hand  the 
depletion  approximation  is  accurate  enouqh  for  most 
practical  purposes.  For  the  potential,  the  difference  in 
results  is  negligible;  for  the  electron  charge 
concentration,  however,  the  exact  n(x)  is  a  smoothly 
varying  function  that  rises  from  10  %  of  Nd  to  90  % 
around  the  depletion  edge  x  *  W,  rather  than  jumping 
discontinuously  from  0  to  Nd  at  x  =  W  for  the  depletion 
approximation.  The  corresponding  potential  change  is  less 
than  2kT/q.  Also  the  electron  concentration  at  x  *  W  is 
about  50  X  of  Nd  in  the  exact  case.  (See  Figure  2.3,  and 
for  further  detail  see  Chapter  3  for  a  one-dimensional 
simulation) 


2.5  PRACTICAL  METAL-SEMICONDUCTOR  DIODE 


The  Schottky  diode  discussed  above  is  an  ideal  one. 
In  this  section  several  practical  departures  from  ideality 
are  discussed. 

2.5.1  Image  Force 

When  an  electron  moves  from  the  semiconductor  to  the 
Schottky  metal,  an  image  force  exists  which  affects  the 
electron  motion.  This  force  lowers  the  potential  barrier 
by  an  amount  (39,40) 


2x  E 
m  max 


(2.23) 


(2.24) 


where  is  the  distance  to  the  peak  of  the  barrier  from 
m 

the  Schottky  contact.  If  the  doping  concentration  is  Nd 
*  1017cm-3,  then  under  zero  bias  E  x  is  of  the  order  of 
106  V/cm,  Xj,  is  of  the  order  of  10  A°r  and  A4>b  is  less 
than  0.05  volt.  This  estimation  suggests  that  the  image 


force  effect  can  be  neglected  in  diode  simulations  if  the 
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2.5.2  Ideality  Factor 

The  current-voltage  relationships  described  by 
equations  2.1  and  2.2  are  for  the  ideal  case.  In 
practice,  no  I-V  characteristics  exactly  obeying  these 
equations  have  been  observed.  An  ideality  factor  n  is 
usually  included  in  the  1-V  relation  to  match  the  theory 
to  experimental  results. 

Jn  =  JQ  exp(-qVQ/kT)  exp (qV/nkT)  [1-exp (-qV/kT) ] 

(2.25) 

where  Jo  is  A*  T2  or  qN„u  E  ,  and  usually  1  <  n  <  1.2. 

an  max 

There  are  several  means  by  which  to  explain  the  factor  n 
added  here,  including  a  voltage  dependent  barrier  height, 
interfacial  layer  effects,  temperature  effects  and  so 
forth.  (For  further  information  see  (39)) 

2.5.3  Tunnelling 

In  reverse  bias,  if  the  potential  barrier  becomes  thin 

enough,  significant  tunnelling  of  electrons  from  the  metal 

to  the  semiconductor  might  take  place.  Also  at  high  doping 

19  ~3 

concentrations,  say  of  the  order  of  10  cm  or  more, 
this  phenomenon  might  be  observed.  In  this  work,  since 
comparatively  low  doping  concentrations  are  used  and  high 
reverse  bias  cases  are  not  of  interest,  the  tunnelling 
effect  is  neglected  in  the  simulations.  For  reference, 
see  (39). 


2.5.4  Interfacial  Layer 

In  the  manufacturing  of  SBDs  it  is  not  easy  to  remove 
the  thin  oxide  layer  often  present  between  the  metal  and 
the  semiconductor.  Such  an  interfacial  layer  lowers  the 
barrier  height  and  gives  some  series  resistance  to  the 
current,  so  that  the  I-V  relationship  departs  from  the 
ideal  case.  Thus  the  ideality  factor  will  change  also. 
Since  the  thickness  of  the  interfacial  layer  might  be  of 
the  order  of  10  A*  ,  to  simplify  the  simulations  in 
two-dimensional  structure,  this  possibility  is  omitted. 

2.5.5  Reverse  Characteristics 

According  to  thermionic-emission  theory  (39),  in 
reverse  bias,  the  current  saturates  at  a  constant  value 

Jsat  =  A*T2  exp(-qVQ/kT)  (2.26) 

However  in  the  real  case,  such  an  ideal  phenomenon  is  not 
seen.  In  fact,  the  current  increases  gradually  with 
voltage.  There  are  several  reasons  which  contribute  to 
this,  including  the  voltage  dependency  of  the  barrier 
height,  the  effect  of  tunnelling,  the  edge  effect,  the 
image  force,  the  carrier  generation  in  the  depletion 
region,  and  so  forth.  Although  the  diffusion  theory 
predicts  an  increasing  reverse  current,  it  needs  some 
modification  to  achieve  a  closer  prediction  to 
experimental  results. 
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When  a  high  reverse  bias  voltage  is  applied,  the 

reverse  current  may  grow  rapidly.  This  sort  of  breakdown 

phonomenon  can  result  from  two  different  mechanisms.  The 

first  one  is  called  Zener  breakdown,  which  operates  at 

lower  voltages;  The  second  one  is  called  avalanche 

breakdown,  which  operates  at  higher  voltages,  beyond  the 

Zener  breakdown.  Zener  breakdown  can  be  thought  of  as  a 

"soft  breakdown"  which  is  due  to  field  ionization  and 

tunnelling.  Usually  such  a  breakdown  phenomenon  occurs  in 

those  diodes  with  heavy  doping  and  under  low  reverse  bias 

voltages.  The  typical  electric  field  is  of  the  order  of 

106  V/cm.  (39,42)  The  avalanche  breakdown  mechanism 

involves  carrier  ionization  and  multiplication  with  high 

reverse  bias  voltages.  (39,42) 

2.5.6  Field  Dependent  Electron  Velocity  And 

Diffision  Coefficient  -  Hot  Electron  Effects 

In  equation  2.13  the  current  is  divided  into  two 

terms.  The  first  is  the  drift  current  qnu  E,  where  u  E  is 

n  n 

the  electron  velocity  Vn  : 

Vn  =  pn E  (2.27) 

The  second  term  is  the  diffusion  current  qDndn/dx. 

In  GaAs,  under  high  fields,  the  electron  velocity,  Vn, 
is  a  nonlinear  function  of  the  field  (43).  According  to 
Thim  (44),  an  analytic  approximation  for  the 


velocity-field  relationship  is 

4  4 

vn  *  IunE  +  vs  (E^}  ]/[1  +  (e^)  1  (2.28) 

o  o 

where  Vs  is  the  saturation  velocity,  Eo  is  the  threshold 
field.  Here  the  temperature  dependent  relation  has  been 
omitted  and  room  temperature  is  assumed.  For  further 
detail  see  (43).  When  the  field  is  low,  equation  2.28 
reduces  to  equation  2.27.  At  very  high  fields,  the 
velocity  saturates  at  a  value  Vs. 

Diffusion  in  GaAs  does  not  obey  the  usual  Einstein 
relation.  However  it  is  possible  to  modify  it  by 
including  an  effective  energy  relaxation  time  in  the 
diffusion  coefficient  (45). 

kT  Vn ^ E  ^  2  2 

Dn(E)  *  -V-  +  f  Te  Vn  <E)  <2*29> 

i 

Here  t  is  the  energy  relaxation  time,  of  the  order 
of  l(T13sec.  for  GaAs.  When  the  field  is  low,  the  second 
term  on  the  right  hand  side  is  negligible  and  equation 
2.29  becomes  the  Einstein  relation. 

For  small  devices  whose  active  region  length  is  less 
than  one  micron,  electron  transport  behavior  is  different 
from  the  discussion  above.  Since  relaxation  time 


processes  for  energy,  charge  and  momentum  are  comparable 
with  the  transit  time,  the  velocity  never  saturates,  but 


rather  increases  with  the  field  (46).  Generally  speaking, 
the  velocity  is  a  function  of  position,  time  and  the 
properties  of  the  external  load,  as  well  as  the  field. 
The  diffusion  coefficient  might  have  a  similar  problem. 
Hot  electron  effects  have  to  be  taken  into  account  in  the 
simulations  presented  here  to  achieve  a  result  closer  to 
reality. 

With  this  as  background,  we  turn  now  to  developing  the 
numerical  techniques  needed  to  solve  the  transport 
equations . 
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COMPUTER-AIDED  DEVICE  SIMULATION 


This  chapter  describes  computer- aided  semiconductor 
device  simulation  techniques  and  presents  a  modified 
decoupled  solution  method  using  a  microcomputer.  The 
one-dimensional  case  is  studied  first.  The  results  are 
evaluated  and  then  the  method  is  extended  to  two- 
dimensional  configurations . 

3.1  PHYSICAL  MODEL 

For  simplicity ,  an  n-doped  majority  carrier  device  is 
assumed,  such  as  a  MESFET,  Schottky  diode,  etc..  The 
characteristics  are  dominated  by  Poisson's  equations  and 
the  electron  current  equation: 


72^  «  SL  (n-N.) 
e  O 


J  /q  *  nv  +  D  Vn 
n  n  n 


(3.1) 


(3.2) 


all  of  the  quantities  were  defined  in  the  previous 

chapter.  For  the  one-dimensional  case  equations  3.1  and 

3.2  are  the  same  as  equations  2.3  and  2.13.  In  this  study 

only  one  and  two-dimensional  behavior  are  considered.  In 

the  two-dimensional  case,  the  parameters  Jn,  vn  and  Vn 

are  vectors  having  components  along  the  x  and  y 
2 

directions;  V  ip  can  also  be  written  as 


The  static  current  constraint 


v-jn  -  o  (3.4) 

is  also  needed  to  solve  the  system.  As  a  matter  of  fact, 

equation  3.2  can  be  combined  with  equation  3.4  to  remove 

the  unknown  current  J  : 

n 

0  -  V« <nVn  +  DnVn)  (3.5) 

At  the  contacts  and  the  edges  of  the  device,  the  boundary 
values  must  be  specified  as  mentioned  in  Chapter  2. 

3.2  NUMERICAL  METHODS 

The  finite  difference*  method  was  the  first  to  be  used 
in  device  simulations.  Five  point  formulas  using 
rectangular  meshes  give  the  famous  centered  difference 
approximation.  In  order  to  solve  the  resultant  system  of 
equations  efficiently,  variable  mesh  sizes  are  considered. 
This  has  been  discussed  by  Adler  (47).  In  many  devices, 
two-  dimensional  behavior  is  restricted  only  to  a  small 
part  of  the  device.  By  separating  the  entire  device  into 
appropriate  portions  depending  on  their  behavior,  and 
defining  suitable  grid  patterns,  much  computer  memory  and 


computation  time  can  be  saved.  In  Adler's  paper,  a  method 
which  allows  terminating  mesh  lines  in  finite  difference 
simulation  was  reported.  This  idea  is  close  to  the  spirit 
of  the  finite  element  method. 

Linear  piecewise  polynomial  functions  applied  on 
triangular  meshes  provide  the  simplest  type  of  finite 
element  formulation.  For  more  accurate  results,  the 
hermite  bicubic  function  method  is  a  better  choice  (24). 
However  this  method  needs  more  boundary  conditions  and 
more  computer  memory  than  the  finite  difference  approach 
and  might  not  be  appropriate  for  small  computer 
simulations.  In  fact,  the  finite  element  method  takes 
more  effort  in  programming  than  does  the  finite  difference 
method.  If  the  geometry  is  not  too  complicated  and  a 
degree  of  approximation  can  be  tolerated,  the  finite 
difference  method  should  be  adguate  for  the  purpose  of 
semiconductor  device  modeling. 

3.2.1  Techniques  Of  Simulation  -  Modified  Decoupled  Method 

As  mentioned  before,  the  finite  difference  method  has 
been  elected  for  simulation.  The  first  step  is  to 
discretize  the  device  by  generating  mesh  lines  over  the 
entire  active  region  under  consideration.  An  example  is 
shown  in  Figure  3.1,  in  which  the  node  numbering  system  is 
also  included.  Mesh  lines  can  be  either  uniform  or 
nonuniform  depending  on  physical  and  geometrical 
considerations.  The  numbering  system  specifies  the 


SCHOTTKY  CONTACT 


boundary  nodes  in  negative  numbers  and  the  nodes  inside 
the  semicoductor  in  positive  numbers.  Values  of  the 
variables  at  the  boundary  nodes,  such  as  those  on  ohmic 
and  Schottky  contacts,  must  be  specified.  At  noncontact 
boundaries,  the  normal  components  of  the  current  and 
electric  field  strength  are  all  equal  to  zero.  However, 
nodes  on  those  surfaces  sure  positively  numbered  and  must 
be  taken  into  account  in  calculation  because  there  are 
unknown  variables  there  as  well  as  at  the  nodes  inside 
the  semiconductor.  Surface  charge  density  can  also  be 
used  as  a  boundary  condition  on  noncontact  boundaries. 

For  each  positive  node  i,  Poisson's  equation  has  a 
corresponding  discrete  version.  For  example  :  in  Figure 
3.1,  if  the  adjacent  nodes  of  i  are  i+l,  i-1,  i+10  and 
i-10  and  all  are  equi -distant  to  i,  then  can  be 
represented  by 

t 

P  (Vi+l+Vi-l-2Vi>  +  ^  <Vi+10+Vi-10-2Vi> 

and  n  simply  by  IT  .  Here  Viand  IT  are  the  space-charge 
potential  and  electron  mobile  charge  concentration  at  node 
i,  respectively.  Therefore,  for  each  node  i,  there  is  a 
discrete  expression  for  Poisson's  equation  with  6  node 
variables:  5  for  potential  and  1  for  charge  concentration. 
This  is  valid  for  most  of  the  nodes  inside  the 
semiconductor .  For  those  nodes  next  to  the  contact. 


however,  one  of  the  node 


the  total  number  of  variables  is  5.  Obviously,  each  node 


on  a  noncontact  boundary  has  fewer  variables  than  those 
inside  the  semiconductor.  Similarly  for  each  node  i,  a 
discrete  current  continuity  equation  (equation  3.5)  can  be 
derived. 

With  suitable  manipulation,  the  entire  family  of 
Poisson's  equations,  after  discretization,  can  be 
expressed  in  matrix  form  as  equation  3.6  and  the  current 
continuity  equation  as  equation  3.7.  Note  that  the  total 
number  of  equations  resulting  from  the  discretization  of 
equation  3.6  and  3.7  is  equal  to  the  total  number  of  nodes 
inside  semiconductor.  For  example,  the  number  is  170  in 
Figure  3.1. 

A  V  =  BN  ♦  C 

P  _  N  =  Q 

(V) 

Here  V  is  a  column  vector  [V1, . ,7^^  r«Presentin9  the 

potential  at  each  node  and  N  =  [N1  . . . N^] t  is  the 

electron  concentration  at  each  node.  The  integer  l 
represents  the  total  number  of  nodes  inside  the 

semiconductor.  A,  B  are  i  x  l  constant  matrices,  C,  Q 

are  constant  l -column  vectors,  P  is  an  l  x  l  matrix  and 
a  function  of  potential.  Note  that  the  current  J  has 
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been  removed  by  equation  3.4.  Now  the  variables  are  V  and 
N.  But  instead  of  solving  for  V  and  N  simultaneously,  in 
equation  3.6,  V  is  treated  as  an  independent  variable,  and 
N  as  a  dependent  variable  with  respect  to  V,  so,  as  a 
matter  of  fact,  equation  3.6  is  solved  in  every  iteration 
loop  with  equation  3.7  as  an  auxiliary  function. 

The  Newton-Raphson  method  (30)  starts  from  a  truncated 
Taylor's  expression.  Concentrating  on  equation  3.6,  after 
some  manipulation  it  becomes: 


(AVi+1-BN1+1-C 


-C)  «  (AVi-BNi-C)  +  A-B  [— ]  AV1  0 

dV 


Vi+1  m  +  AV^ 


(3.8) 


(3.9) 


where  i  =  1,  2,  3,  .  represents  the  index  of 

dN 

iterations.  The  quantity  [-3]  is  an  l  x  l  tangential 

dV 

matrix  or  Jacobian  matrix’ (30)  : 


3 


AV  can  be  evaluated  from  equation  3.8  with  the  result 


AV1 


s  ii-1 


-  (A 


"B[^J  | 

dV  ) 


(AVi-BN1-C) 


(3.10) 


In  the  above  equation  two  terms  are  still  unknown,  that  is 


dN  1 

Frl 


and 


N1 


dV 


These  values  can  be  found  from  equation  3.7: 


n  =  [p  _  r1  q 


(3.11) 


CV) 


and  also  by  differentiating  it  with  respect  to  V, 


[-=  (P  N)  ]  **  0 
dV  (V) 


(3.12) 


The  left  hand  side  of  equation  3.12  can  be  separated  into 
two  terms. 


_  rdN.  .  _  l 

P  -  [— ]  +  (Yik) 

(V)  dV  j ,k*l 


(3.13) 


t  O  - 


where 


jk 


him# 


v.  and  n. 


m 
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are  elements  of  P 


(V) 


,  V  and  N  respectively. 


3.13  may  be  rearranged  to  read: 


Equation 


dV 


-  [P  _  3 

(v) 


-1 


[Y-jk] 

3K  j /k=l 


(3.14) 


This  completes  the  evaluation  of  AV1  for  equation  3.10. 
From  equation  3.9,  V1+1  can  be  calculated  and  equations 

3.9  and  3.10  are  ready  for  next  iteration  loop  i+1 . 

To  start  the  calculation,  the  initial  V  value  is 

arbitrarily  given,  say  V  *  V°  «  [0, . ,0]fc.  From 

equation  3.7,  N  can  be  obtained,  say  N?  which  is  used  in 
the  right  hand  side  of  equation  3.6.  Then  equation  3.6 
can  be  solved  directly  in  terms  of  V.  This  V  is  used  as  V1 
and  the  aforementioned  iteration  procedure  starts . 

Usually  it  takes  3  to  6  iterations  to  obtain  a  converged 
solution  for  V.  The  solution  of  N  follows  through 
equation  3.11. 

The  numerical  approach  described  above  can  be 
summarized  as  follows  :  Eliminate  N  in  equation  3.6  by 
equation  3.11  and  obtain 


AV  *  B[P  ]  -1  Q  +  C  (3.15) 
(V) 

This  is  a  system  of  equations  with  variable  V  only. 
Theoretically  it  is  possible  to  solve  such  nonlinear 
system  using  the  Newton-Raphson  method  directly.  However, 


difficulties  arise  from  manipulation  of  the  inverse  term 
on  the  right  hand  side  of  equation  3.15.  This  explains 
why  it  takes  so  many  steps  to  process. 

From  the  description  above  it  is  evident  that  the 
largest  memory  space  occupied  by  one  term  is  an  l  x  l 
matrix  storage.  Since  there  is  only  one  such  matrix 
needed  for  processing  at  any  one  time,  all  of  those  terms 
not  in  use  can  be  stored  on  disk.  For  manipulation  of 
equation  3.14,  there  are  two  l  x  I  matrix  multiplications. 
Since  each  matrix  may  occupy  more  than  85  %  of  total  CPU 
memory  space,  one  of  them  should  be  stored  on  the  disc, 
and  the  matrix  operation  done  column  by  column. 

3.3  COMMENTS 

In  this  section  some  practical  considerations  about 
this  simulation  are  disscussed. 

3.3.1  Boundary  Conditions 

It  is  necessary  to  consider  the  boundary  conditions 
when  the  simulation  is  carried  out.  At  ohmic  contacts, 
usually  the  ideal  case  is  assumed,  i.e.  the  carrier 
concentration  is  equal  to  the  doping  concentration  of  the 
adjacent  semiconductor.  At  the  Schottky  contacts,  in  the 
forward  bias  case  the  carrier  concentration  is  assumed  to 
be  given  by  equation  2.16.  In  reverse  bias,  at  the 
Schottky  contacts,  since  the  effective  built-in  potential 
becomes  VQ  +  Vr  the  following  expression  is  more  suitable 


Here,  Vr  is  the  magnitude  of  the  reverse  bias  voltage. 
Also  at  the  edges  of  the  device,  it  is  assumed  that  the 
first  order  partial  derivatives  of  both  carrier 
concentration  and  space-charge  potential  in  the  normal 
direction  are  zero. 

3.3.2  Current  Continuity  Equation 

For  static  simulations,  the  unknown  current  is 
eliminated  by  combining  equations  3.2  and  3.4  to  obtain 
equation  3.5.  This  makes  the  equation  a  second  order 
partial  differential  equation  as  well  as  a  Poisson's 
equation.  So,  in  fact,  two  sets  of  boundary  values  are 
required  for  simulations. 

If  the  maximum  field  everywhere  in  the  device  is  low, 
the  Einstein  relationship  holds,  and  the  simulation  is 
straightforward.  If,  however,  the  field  is  high  enough  to 
drive  carriers  into  the  nonlinear  parts  of  the 
velocity- field  curve  and  the  diffusion-field  curve  for  a 
significant  part  of  the  carrier's  trajectory,  (in  case  of 
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1.  Equations  2.16  and  3.16  are  still  used,  due  to  lack 
of  better  boundary  conditions  available. 

2.  Solve  Poisson's  equation  and  the  continuity  equation 

numerically,  in  one  dimension,  taking  hot  electron 
effects  and  the  value  of  the  current,  i.e.  equation 
2.2,  into  account,  to  obtain  the  carrier 

concentration  at  the  Schottky  contact  for  each 
applied  voltage. 

The  first  choice  is  simple  but,  in  two-dimensional 
simulations,  the  total  current  calculated  would  be 
unreasonable  under  certain  conditions.  However,  in 
two-dimensional  simulation  of  MESFETs  for  example,  usually 
the  metal  gate  is  reverse  biased,  and  negligible  current 
is  carried  by  the  gate.  So  equation  3.16  and  zero  current 
density  are  both  assumed  without  loss  of  generality  (20). 
The  second  choice  might  be  good  for  fine  meshes.  However, 
the  memory  capacity  of  the  small  computer  used  here  limits 
the  smallest  mesh  size,  and  so  in  some  cases  unreasonable 
carrier  concentrations  at  the  Schottky  contact  are 
obtained. 

3.3.3  Program  And  Flow  Chart 

Since  one  wants  to  reserve  most  of  the  computer  memory 
for  matrix  operations,  the  program  is  separated  into 
several  subprograms .  All  the  data  not  in  use  are  stored 
on  disk.  Some  matrix  variables  and  array  variables  can 
represent  different  items  at  different  times  to  minimize 


the  number  of  variables  required.  Caution  must  be 
exercised  in  programming  equation  3.10,  especially  in 
calculating  the  Jacobian  matrix  as  in  equation  3.14,  to 
prevent  some  error  in  differentiating  each  element. 

Information  regarding  each  node  or  element  is  stored 
on  the  disk,  including  the  adjacent  node  or  element 
number,  the  distance  between  each  of  the  nodes  and  some 
identification  to  distinguish  between  nodes  inside  the 
active  region  and  those  on  the  boundary.  The  boundary 
values  can  be  adjusted  automatically  in  the  main  program 
for  each  applied  bias.  So  the  program  can  be  used  to 
simulate  different  structures  if  the  input  information  is 
properly  assigned. 

To  start  the  simulation,  the  value  of  the  space  charge 
potential  at  each  node  is  assumed  under  the  zero  bias 
condition.  Since  the  current  is  zero  in  this  case,  one 
can  solve  for  the  carrier  concentration  from  current 
continuity  and  carry  out  an  iteration  procedure  until  it 
converges.  Usually  for  a  simulation  containing  80  nodes, 
three  to  six  iteration  loops  are  common  and  it  takes 
approximately  1.5  to  2.5  hours  to  obtain  the  results  for 
one  bias  point  on  the  device  I-V  curve,  using  an  HP  9825B. 

The  flow  chart  for  the  process  is  shown  in  Figure  3.2 
and  a  sample  program  is  given  in  the  appendix.  Figure  3.3 
shows  the  steps  needed  to  operate  the  program.  An 
illustration  of  the  device  configuration  with  mesh  lines, 
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and  the  nodes  numbered  is  shown  in  Figure  3.1.  To  store 
node  information,  the  first  program  in  the  appendix  can  be 
used  as  an  example.  Boundary  values  can  be  defined  in  the 
main  program  such  that  after  the  information  for  each  node 
is  read  in,  its  boundary  value  will  be  given  automatically 
before  the  calculation  proceeds,  as  shown  in  lines  17-28 
of  Program  3  in  the  appendix.  The  final  results  of  V  and 
N  for  each  applied  bias  should  be  stored  on  disk  and  also 
printed  out  in  the  same  form  as  node  numbering 
configuration.  A  typical  program  for  this  is  shown  by 
Program  7  in  the  appendix.  Finally  the  I-V  curve  can  be 
plotted.  In  the  next  section,  considerations  needed  to 
set  mesh  sizes  are  discussed. 

3.3.4  Mesh  Spacing 

In  semiconductor  device  simulations,  it  is 
recommended  to  restrict  the  mesh  spacing,  h,  to  be  no  more 
than  the  Debye  length: 

w  -  (^r 

i.e. 

2  2  1/2 

Lpebye  ^  [(dx)  +  {dy)  ]  ”  h  (3*18) 

where  dx  and  dy  are  the  mesh  spacing  components  in  the  x 
and  y  directions,  respectively.  If  this  limitation  is 


ignored,  the  solutions  may  represent  physically 
meaningless  results  (24).  This  restriction  is  mainly  for 
those  cases  where  the  Einstein  relationship  holds  or  at 
least  for  deviations  not  far  from  it.  If  some  nonlinear 
velocity- field  and/or  diffusion-field  curves  are  used,  the 
maximum  spacing  might  be  modified  as  (24) 


Vi  Ay  <  2D 


(3.19) 


where  V  is  the  saturation  velocity,  and  av  is  the 
sat 

potential  change  between  two  adjacent  nodes.  Also,  from 
the  discussion  above,  since  the  maximum  spacing  is 
limited,  the  maximum  device  area  which  we  can  model  is 
also  limited  because  of  the  computer  memory  size. 

As  an  example  to  obtain  some  understanding  about  the 
minimum  spacing,  consider  the  conventional  SBD  of  Figure 
2.3.  Under  zero  bias,  from  equation  2.18,  Poisson's 
equation  becomes 


0  -  |Nd|exp  (*-V0)J-  lj 


(3.20) 


At  any  point  x  located  in  the  area  between  the  depletion 
edge  W  and  the  end  of  the  epi-layer  L,  is  very  close  to 
Vo.  Using  Taylor's  expansion  in  equation  3.20  yields 


(3.21) 


Let  c  s  ^  -  VQ,  so  equation  3.20  cam  be  rewritten  as 

7§  *  I  el «  <3-22> 

ax 

If  the  Newton-Raphson  method  is  used,  the  truncation  error 
is 

^h2sf| 

A  dx^ 

Assuming  the  mesh  spacing  h  is  no  less  than  the  Debye 
length,  then 

1  h2  -  I  LDebye  I  Nd  3#  5  (3’23) 

-  5 

In  other  words,  the  error  could  be  larger  than  100  %.  For 
two  or  three-  dimensional  structures,  similar  results  can 
be  expected . 

The  argument  given  above  for  minimum  spacing  is  no 
longer  valid  inside  the  depletion  region,  since  in  this 
case  equation  3.20  becomes 


instead  of  equation  3.21.  It  can  be  shown  that  the 
truncation  error  is  -kT/q,  which  is  comparatively  small 
compared  with  the  potential  value  V  corresponding  to  it. 
Therefore ,  wider  spacing  is  allowed  inside  the  depletion 
region  than  outside  it.  In  general,  the  Debye  length  is  a 
good  measure  for  the  mesh  spacing.  Caution  should  be 
exercised  when  the  mesh  spacing  is  chosen  to  be  larger 
than  Debye  length. 

If  the  Einstein  relationship  is  not  applied,  then 
suitable  mesh  spacing  related  to  the  current  equation 
should  also  be  considered.  Equation  3.19  is,  in  general, 
a  basic  criterion  for  that.  Generally,  if  the  mesh 
spacing  is  taken  to  satisfy  both  equations  3.18  and  3.19, 
reasonable  results  should  be  expected.  From  simulation 
cost  considerations,  it  is  possible  to  relax  the  spacing 
criterion  in  the  interior  a  certain  amount,  but  keeping  a 
finer  mesh  spacing  on  the 'boundaries. 

3.4  ONE-DIMENSIONAL  SIMULATION 

Since  the  conventional  Schottky  diode  is  reasonably 
well  understood,  it  can  be  used  for  comparison  with  the 
one-dimensional  simulation  results.  If  they  are  similar, 
then  some  confidence  can  be  placed  on  the  results  of  the 
two-dimensional  simulations.  In  this  section  we  will 
consider  only  one-dimensional  cases  and  examine  the 
details  of  the  static  solution  in  a  conventional  GaAs  SBD. 
The  results  shown  in  Figure  3.4  show  the  effect  of  varying 
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Figure  3.4  Typical  results  for  a  conventional  Schottky- 
barrier  diode,  (a)  potential  (b)  carrier 

concentration  with  ■  10*^  ca"\  barrier 

height  ■  0.9  V,  un  ■  6875  an^V"*sec"*,  Vg  » 

1.1  x  lO^cm/sec,  E  ■  3.2  KV/aa,  Starting  from 
the  solid  line  (Einstein  relation) ■  0,  1, 

3,  5,  7,  9  x  10‘^sec,  depletion  length  ■  0.11 


the  energy  relaxation  time.  From  this  graph  and  the 
results,  one  observes  that: 


1.  For  a  fixed  bias  voltage,  the  field  strength  at  the 
Schottky  contact  is  approximately  a  constant,  no 
matter  what  kind  of  field  dependent  velocity  and 
diffusion  relationships  are  taken.  Recall  that 


This  result  can  be  understood  as  follows:  From 
equation  2.3, 

(i£)L  •  (ix)0  ■  J  2  I  (n_Nd)  §&  (3-25) 

*"  o 

Assuming  the  diode  length  L  is  greater  than  the 
depletion  length  W,  then 


(3.26) 


Also,  n  is  negligible  inside  the  depletion  region. 
Bence 


li/2 

Nd<Vv)J 


(3.27) 


where  the  right  side  of  the  equation  above  is  a 
constant  for  any  given  V. 


2.  For  the  potential  curves,  near  the  Schottky  contact, 

all  the  curves  are  close  together;  near  the  edge  of 
depletion  region  they  deviate  from  one  another  by  no 
more  than  2kT/q.  Also  they  are  close  to  the 
potential  curve  found  from  the  depletion 

approximation . 

3.  For  electron  concentration,  all  the  curves  meet  near 

the  edge  of  depletion  region  with  electron 

concentration  approximately  eqral  to  50  %  of  Nd  . 

The  total  charge  Q  in  the  depletion  region  for  each 

curve,  is  equal  to  MjiJ)  multiplied  by  the  diode 

ax  0 

area,  which  is  the  same  as  in  the  depletion 
approximation. 

Thus,  although  the  detailed  shape  of  the  depletion 
region  differs  from  that  predicted  by  the  depletion 
approximation,  so  little  charge  is  displaced  that  only 
negligible  effects  on  the* field  (and  hence  the  potential) 
are  observed.  This  displaced  charge,  however,  will  be 
seen  to  have  more  significance  when  we  wish  to  explain  the 
properties  of  the  Gap  diode  in  the  next  chapter. 


3 . 5  TWO-DIMENS IONAL  S IMULATION 
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With  64  K-bytes  of  memory,  the  maximum  number  of  nodes 
available  for  calculations  is  about  81.  Under  this 
limitation  one  has  to  suitably  define  the  meshes  over  the 
entire  area.  All  the  node  information  is  stored  on  the 
disk  along  with  subprograms,  intermediate  data,  etc.. 
Most  of  the  computing  time  is  used  in  calculating  the 
inverse  matrices  of  equations  3.10  and  3.14. 

If  the  relationship  between  Dn  and  \in  is  close  to  the 
Einstein  relationship  or  Dn  /un  -  constant,  the  following 
arrangement  can  be  applied  to  simplify  the  simulation:  The 
current  equation  may  be  rewritten  as: 

V1  jj 

Jn/q  “  Dn  exP(D^  Vln  exP<“  5^  *>]  (3.28) 
n  n 


If 

n  Wn 

if-  exp  (-  R—  =4* 
d  n 

is  defined,  equations  3.1  and  3.28  become: 


V2*  -  |  Na 


(<{>  exp  ip)  -  1] 
n 


(3.29) 


and 


Jn/q  “  Dn  exp  (^  ip)  Nd 

n 


V<(> 


(3.30) 
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respectively.  Now  the  varieties  are  <t>  and  \p  instead  of  n 
and  \p .  The  simulation  procedures  are  the  same  as  that 
described  previously.  The  advantage  of  this  is  that  it  is 
easier  to  handle  the  boundary  values  this  way  and  the 
calculation  is  simplified.  In  most  of  the  cases,  it  is 
possible  to  adjust  the  value  of  $  to  be  in  the  range  of  0 
to  1,  in  order  to  obtain  accuracy  and  convergence.  Also, 
the  current  has  the  same  direction  as  v<p  >  so  one  can 
easily  predict  the  current  direction  at  each  node.  This 
also  helps  to  check  if  the  simulation  results  are 
physical . 

In  the  next  two  chapters  two  simulation  examples  and 
their  results  are  shown. 
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4.  GAP  CONTROLLED  METAL- SEMI CONDUCTOR  DIODE 


In  this  chapter,  a  new  metal-semiconductor  device  is 
described  whose  I-V  characteristic  is  controlled  by  a 
geometric  gap  between  adjacent  Schottky  metal  regions. 

4.1  CONFIGURATION 

Consider  the  cross-sectional  view  shown  in  Figure  4.1. 
The  device  represented  there  consists  of  an  n-type 
epitaxial  layer  of  GaAs  on  an  N+  conducting  substrate. 
Stripes  of  a  Schottky  metal  are  deposited  either  on  the 
surface  or  in  notches  in  the  epi-layer.  Thus,  a  gap 
exists  between  adjacent  Schottky  barrier  regions  exposing 
the  n-layer.  Finally,  an  ohmic  contact  is  established 
over  the  entire  surface,  connecting  the  Schottky  barrier 
regions  and  the  n-layer  filled  gap  in  parallel.  Note  that 
this  cell  can  be  repeated  across  the  surface  of  a  wafer. 
We  term  this  configuration  the  Gap  diode. 

4.2  COMPUTER  SIMULATIONS  AND  ANALYSIS 


The  parameters  used 
listed  as  follows: 
built-in  potential 
doping  concentration 
energy  relaxation  time 
saturated  Velocity 
critical  field 
diode  height 


in  the  device  simulation  are 

<f>b  =  0.72  V 

N,  =  2xl015  cm-3 
a 

Te  =  10"13  sec 
V  =  2xl07  cm/sec 

S 

Eq  =  3.2  KV/cm 
12  =  0.81  um 
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Ohmic  Contact 

//////////////, 

Gap 

§|§||| 

n 

N+-  Substrate 

Figure  4.1  Cross-section  of  the  Gap  Diode.  Typical  device 

;  gap  width  -  1  ym 
diode  area  ■ 


2  x  1015c»"3 


parameters: 
epi-layer  thickness  -  1  ym 
8  x  10-6cm"2. 


gap  width 

di 

=  0.72,  1.08, 

aspect  ratio 

V1! 

=  8/11 

Schottky  height 

d2 

=  0.36  iim 

diode  area 

Vs 

=  8xl0-6  cm2 

depletion  length 

Mo 

=  0.706  pm 

The  configuration  is  shown  in  Figure  4.2. 

The  simulation  procedures  and  techniques  have  been 
described  in  Chapter  3.  Figure  4.3  shows  the  mesh  lines 
and  the  Schottky  contacts  denoted  by  filled  dots.  Note 
that  the  entire  area  can  be  separated  into  two  symmetrical 
parts.  Therefore,  only  the  upper  half  is  simulated.  The 
open  circles  represent  the  ohmic  contacts  and  the  crosses, 
the  noncontact  boundaries.  This  information  is  entered  in 
Program  1  of  the  appendix.  For  example:  In  line  41,  "I" 
denotes  nodes  inside  the  semiconductor.  In  line  40,  "B" 
denotes  nodes  on  the  boundaries.  Following  this,  line  61 
asks  whether  the  node  is • on  the  Schottky  contact  or  the 
ohmic  contact,  etc.  In  the  same  program,  the  following 
information  is  also  given:  the  mesh  spacing  (lines  25-28, 
57-60),  the  adjacent  nodes  (lines  51-54),  geometry 
structure  control  (lines  35,  90-111),  and  material 
parameters  including  doping  concentration,  barrier  height 
and  electron  mobility  (all  in  line  32).  Finally  all  node 
information  is  stored  in  a  file  named  Cnode  (lines  71-74), 
and  all  others  in  Cinp  (line  112)  and  CJ#  (line  127),  here 
file  CJ#  is  for  current  density  input  and  output  control. 


Figure  4.3  Mesh  lines  used  in  simulation. 

Only  the  upper  half  including  the  symmetry 
line  is  involved  in  the  simulation. 
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Boundary  value  assignment  is  set  in  the  main  program 
(lines  17-26,  Program  3  in  the  appendix). 

A  typical  printout  of  the  computer  results  is  shown  in 
Figure  4.4.  The  first  part  is  the  potential  value  at  each 
node.  The  second  is  the  coreesponding  carrier 
concentration  normalized  to  doping  concentration.  The 
printout  has  been  designed  to  have  the  same  appearance  as 
Figure  4.3.  Also,  the  current  density  at  each  node  on 
both  sides  is  given  to  help  in  understanding  the  current 
conduction.  Finally,  the  average  current  density  on  the 
bottom  line  multiplied  by  the  diode  area  is  the  total 
current.  The  average  current  density  values  on  the  input 
and  output  sides  also  serve  as  a  check  of  the  quality  of 
the  simulation.  Notice  that  the  round-off  error  might 
limit  the  precision. 

From  these  numerical  results,  iso-electron 
concentration  and  equi -potential  plots  can  be  drawn,  as 
shown  in  Figures  4.5  and  4.6.  Note  that  the  curves  near 
the  corners  of  gap  area  adjacent  to  the  Schottky  metal  are 
not  precise,  because  of  the  uncertainty  for  the  boundary 
values.  These  two  figures  give  some  insight  into  how  the 
diode  works. 

The  principal  of  operation  is  as  follows.  At  zero 
bias  the  Schottky  metals  deplete  the  semiconductor  both 
below  the  contacts  themselves  and  in  the  gap  between  them. 
By  adjusting  the  width  of  the  gap  and  the  depth  of  the 
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Figure  4.4  Typical  printout  of  computer  results 


Figure  4.5  Normalized  electron  concentration  contours 
for  one-depletion-length  Gap  diode. 

(a)  forward  bias  :  0.3  V  (b)  zero  bias 
(c)  reverse  bias  :  -0,3  V. 


notch,  the  turn-on  voltage  for  a  given  doping  and  choice 
of  metal  can  be  set.  For  example.  Figure  4.5  (b)  shows 
contours  of  equal  free  charge  concentration,  n,  normalized 
to  the  doping  concentration,  N^,  with  no  bias  applied. 
Note  the  heavily  depleted  regions  under  the  Schottky 
contacts  and  up  into  the  gap  as  well.  In  reverse  bias. 
Figure  4.5  (c),  the  depletion  regions  punch  through  to  the 
cathode  contact  and  the  neck  between  the  depletion  regions 
moves  down  from  the  gap.  Reverse  leakage  current  is 
supported  by  the  large  number  of  free  carriers  in  the  gap 
and  in  the  neck  between  the  depleted  regions.  In  forward 
bias.  Figure  4.5  (a),  the  depletion  regions  are  retracted 
toward  the  metal  contacts,  thus  exposing  an  ohmic  channel 
for  charge  carrier  conduction.  Notice  that  this  occurs 
for  forward  bias  well  below  that  needed  to  draw  current 
through  the  Schottky  contacts.  This  cam  be  seen  in  the 
current  density  results  in  Figure  4.4,  which  also  shows 
the  effect  of  current  spreading  at  the  cathode  contact. 

It  can  be  understood  through  the  one-dimensional 
simulation  in  Chapter  3  that  the  edge  of  the  depletion 
region  is  not  so  sharp  as  predicted  by  the  depletion 
approximation.  The  displaced  charge  between  the  two 
depletion  regions  in  the  gap  explains  why  forward 
conduction  through  the  gap  is  much  easier  than  through  the 
Schottky  contact.  This  also  explains  why  the  reverse 
current  is  much  higher  than  in  a  conventional  Schottky 


I 


contact.  However ,  the  reverse  current,  in  general,  is 
nuch  less  than  the  forward  current,  since  the  depletion 
region  spreads  into  the  gap  in  reverse  bias. 

Current  -  voltage  characteristics  for  the  Gap  diode 
are  shown  in  Figure  4.7.  Here  the  parameter  is  the  gap 
width  normalized  to  the  zero  bias  depletion  length.  Wo. 
The  right-most  curve  is  the  result  for  a  conventional  SBD. 
Note  that,  as  expected,  the  turn-on  voltage  can  be  set 
from  nearly  the  built-in  potential  down  to  nearly  zero. 
However,  as  the  gap  is  widened  the  reverse  current  is 
increased. 

Figure  4.8  shows  the  corresponding  I-V  curves  in  log 
scale.  For  comparison,  a  pure  resistor  which  is 
equivalent  to  a  gap  diode  with  a  very  wide  gap  length,  say 
four  depletion  lengths  or  more,  is  presented.  Note  that 
the  Schottky  diode  Z-V  curve  is  a  straight  line  in  log 
scale.  As  the  gap  length* increases,  the  curvature  changes 
correspondingly.  When  the  gap  width  exceeds  two  depletion 
lengths  the  I-V  curve  will  converge  to  that  of  a  resistor. 
Thus,  the  shape  of  the  Gap  diode  I-V  curve  is  also 
controllable  by  adjusting  the  gap  length.  For  different 
applications,  different  requirements  on  the  tum-on 
voltage  and  I-V  curvature  might  be  preferred.  Then  the 
optimization  of  the  gap  length,  designed  to  satisfy  a 
special  application,  will  deserve  to  be  studied  in  further 
detail . 
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Figure  4.9  shows  predicted  I-V  curves  using  the  depth 
of  the  notch  d2  as  the  control  variable.  Thinner  notches 
correspond  to  lower  tum-on  voltages.  Note  that  even  for 
d  =0,  the  I-V  curve  is  still  controllable  by  adjusting 
the  gap  length.  This  is  the  planar-type  Gap  diode. 

The  capacitance  can  be  calculated  from 


c  = 


dQ 

1 

dV 

AV 

I  T,  (Nd-Ni) 
1  =  1 


V+AV 

f 

V 


111213 


(4.1) 


Here,  l  and  Ni  are  defined  in  Chapter  3.  l^l^  is  the 
volume  of  the  bulk  semiconductor  region.  The  predicted 
capacitance  -  voltage  relation  is  shown  in  Figure  4.10. 
Before  tum-on,  the  capacitance  of  the  Gap  diode  is  lower 
than  that  of  the  same  size  Schottky  diode.  This  is 
because  most  of  the  current  goes  through  the  gap  instead 
of  the  Schottky  metal  contact.  When  the  bias  voltage 
approaches  the  tum-on  voltage,  which  is  approximately 
0.46  volts  in  this  case,  the  capacitance  increases 
rapidly.  For  the  Gap  diode,  the  current  conduction  can  be 
considered  as  being  controlled  by  an  equivalent  potential 
barrier  in  the  gap.  Its  value  is  close  to  the  tum-on 
voltage.  Recall  that  the  capacitance  is  defined  in 
equation  2.12.  When  the  applied  bias  approaches  the 
tum-on  voltage,  W  becomes  very  small  and  the  capacitance 
increases  very  quickly.  Note  that  in  Figure  4.10  the  SBD 
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Figure  4.10  C-V  curve  of  a  Gap  diode  with  a  one  depletion  length 
gap  in  comparison  with  a  corresponding  Schottky- 
barrier  diode. 
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was  assumed  to  have  infinite  length.  However,  the  length 
of  the  Gap  diode  is  approximately  one-depletion  length. 
Therefore,  in  the  reverse  bias  case,  the  C-V  curve  of  the 
Gap  diode  approaches  a  low  value  rapidly  due  to 
punch- through  and  the  conduction  modulation  effect  of  the 


4.3  EXPERIMENT 

Experiments  with  this  new  device  have  been  reported  by 
Teng  et  al.  (37).  Figure  4.11  shows  the  experimental  I-V 
characteristics  for  diodes  with  seven  different  gap 
widths.  The  right  hand  curve  is  for  a  Schottky  diode  with 
a  30  iim  diameter,  the  same  as  that  of  the  Gap  diodes. 
Note  the  ranking  similarity  to  the  predicted  I-V  curves 
of  Figure  4.7. 

4.4  COMMENTS 

1.  Effort  has  been  expended  to  find  an  analytical 
expression  for  the  I-V  relation.  One  way  to  treat  this 
problem  is  to  apply  the  potential  barrier  concept  on  the 
gap  region  and  consider  the  entire  diode  as  a  parallel 
combination  of  the  Schottky-barrier  region,  rectifing  gap 
region  and  perhaps  the  ohmic  gap  region  if  the  gap  is 
comparatively  large.  There  is  no  difficulty  in  finding 
the  I-V  relation  for  the  Schottky  barrier  region. 
However,  because  of  the  narrow  gap  nature,  it  is  difficult 
to  obtain  an  suitable  expression  for  the  current  through 
the  gap. 


Figure  4.11  Experimental  I-V  characteristics  for  diodes  with 
seven  different  gap  widths.  Right  hand  I-V 
curve  is  for  pure  Schottky-barrier  diode. 


2.  In  the  gap,  since  both  Schottky  contacts  have  the 
same  potential,  therefore,  right  underneath  the  ohmic 
contact  there  should  be  some  depletion  area  on  both  sides. 
Because  of  the  properties  of  Poisson's  equation,  the 
transitional  depletion  edge  may  stretch  out  to  as  much  as 
two  Debye  lengths.  Also,  because  of  the  interaction 
between  the  two  adjacent  depletion  regions,  the  total 
controllable  gap  length  might  be  as  much  as 


W  +  2L_  .  „  +  6  +  2Ln  .  +  V?  *  3  V* 

o  Debye  Debye  o  o 


(4.2) 


where  6  is  of  the  order  of  one  Debye  length.  Note  that 


Wo/LDebye 


.r^l 

LkT/QJ 


*  5.27 


(4.3) 


if  Vo  =  0.72  V. 

This  is  close  to  what  is  observed  experimentally. 

3.  We  have  demonstrated  a  new  device  principle:  It  is 
possible  to  produce  depletion  of  carriers  under  an  ohmic 
contact  region.  By  this  means,  a  large  range  of  turn-on 
voltages,  from  that  of  the  pure  Schottky  to  nearly  zero 
volts,  can  be  obtained.  This  observation  can  be  applied 
to  other  two-  or  three-  terminal  structures  which  may 
yield  microwave  and  mm-wave  devices  with  new  and  desirable 
properties . 
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5-  GEOMETRICALLY  CONTROLLED  ANTISYMMETRIC 
METAL-SEMICONDUCTOR  DIODE 


In  this  chapter  the  design  of  an  antisymmetric  diode 
is  reported.  The  computer  results  show  that  the  shape  of 
its  I-V  curve  can  be  adjusted  as  well  as  its  tum-on 
voltage.  This  new  device  has  possible  high  frequency 
applications . 

In  1974  Schneider  and  Snell  (48)  and  Cohn,  Degenford 
and  Newman  (38)  reported  simultaneously  a  new  mixer,  named 
the  antiparallel  mixer  diode  pair.  In  this  device, 
because  of  its  antisymmetric  configuration  (See  Figure  5.1 
(b)),  the  I-V  curve  is  symmetrical  with  respect  to  the 
origin.  If  a  local  oscillator  sinusoidal  voltage  is 
applied  to  this  device,  there  will  be  two  switching  cycles 
for  each  cycle  of  the  local  oscillator  voltage.  This  is 
double  that  of  the  single  diode  mixer  case.  It  has  been 
shown  (38)  that  the  total  current  of  the  antiparallel 
diode  pair  contains  only  frequencies  mf  +  nf  for  which 
m  +  n  is  an  odd  integer;  f  and  f  are  the  frequencies  of 
the  local  oscillator  and  the  signal,  respectively.  On  the 
other  hand,  the  single  diode  mixer  results  in  a  diode 
current  having  all  frequency  components.  This 
characteristic  property  suggests  that  such  a  configuration 
has  the  advantage  as  a  subharmonic  mixer,  because  it  can 
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be  pumped  at  half  of  the  frequency  required  for  the  local 
oscillator  in  a  single  diode  mixer. 

The  natural  suppression  of  the  fundamental  and  other 
mixing  products  is  lessened  due  to  the  slightly 
nonidentical  electrical  properties  of  the  two  diodes  for 
the  antiparallel  diode  mixer  pair,  such  as  differences  in 
saturation  currents  and  I-V  curvatures  (38).  Furthermore, 
although  the  diodes  may  be  identical,  since  they  are 
physically  separated  their  circuit  environments  may  not  be 
identical.  Therefore,  the  overall  electrical 
characteristics  of  the  antiparallel  mixer  diode  pair  may 
not  be  symmetric.  This  also  results  in  the  incomplete 
suppression  of  undesired  mixing  terms.  So  the  conclusion 
is  that  one  needs  good  electrical  symmetry  as  well  as 
mechanical  symmetry  for  proper  device  operation  (49).  Due 
to  practical  considerations,  this  is  especially  hard  to 
accomplish  at  mm-wave lengths  (50). 

Also  mentioned  by  Cohn  et  al.  (38),  the  suppressed 
current  components  circulate  within  the  loop  formed  by  the 
two  diodes.  According  to  Kerr  (49),  the  magnitude  of  the 
loop  inductance ,  as  seen  by  currents  circulating  through 
the  two  diodes,  strongly  affects  the  subharmonically 
pumped  mixer  per f ormance ,  such  as  conversion  loss  and 
noise  figure  in  mm-wave  applications.  Therefore, 
adjusting  this  loop  inductance  to  optimize  the  performance 
is  an  important  topic  to  be  considered. 
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In  1978,  Cardiasmenos  reported  a  successful  technique 
to  fabricate  such  a  device  with  a  5-dB  DSB  noise  figure  at 
94  GHz  ( 50 ) .  These  planar  type  devices  are  used  in 
suspended  stripline  substrates.  Microprocessor-controlled 
machinery  replaces  the  traditional  trial  and  error 
fabricating  and  trimming  work  and  achieves  highly 
reproducible  results.  The  entire  fabrication  procedure  is 
very  complicated  and  this  is  believed  to  be  the  first 
successful  report  on  the  mass  production  of  the 
antiparallel  mixer  diode  pair. 

In  order  to  explore  alternative  means  of  realizing  the 
desired  I-V  characteristics,  we  consider  an  innovative 
design  for  such  device  in  the  following: 

5.1  CONFIGURATION  AND  PHYSICAL  ANALYSIS 

Consider  the  geometry  of  Figure  5.1  (a).  A  Schottky 
metal,  such  as  A.  is  deposited  over  a  portion  of  an  N+ 
substrate,  and  N-doped  GaAs  layer  is  grown  around  and  over 
it.  After  the  desired  diode  thickness  is  obtained  a 
second  Schottky  metal  is  deposited  over  the  antisymmetric 
side  and  finally  an  ohmic  contact  is  applied  to  the  entire 
upper  surface,  completing  the  device.  We  call  this  device 
the  antisymmetric  diode.  The  equivalent  circuit  of  the 
device  is  shown  in  Figure  5.1  (b)  which  is  the 
antiparallel  diode  pair.  The  technology  needed  to 
fabricate  such  a  device  is  the  same  as  that  used  in  the 
permeable  base  transistor  (51). 
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(a)  Antisymnetric  diode  configuration. 


(b)  Equivalent  circuit  -  antiparallel  (back-to-back) 

diode  pair. 


Figure  5.1  Antisymmetric  diode  and  its  equivalent  circuit. 


Such  a  structure  should  be  highly  symmetric  and 
uniform  so  that  well  balanced  I-V  curves  should  be 
obtained.  Furthermore,  the  millimeter  wave  parasitic 
circuit  elements  are  all  external  to  the  new  diode  so  that 
circuit  balance  should  be  improved  and  impedance  matching 
simplified.  In  addition  to  the  mechanical  and  electrical 
advantages  of  this  structure  it  can  be  shown  that,  with 
the  barrier  height  unchanged,  the  tum-on  voltage  can  be 
lowered  to  a  suitable  value. 

5.2  COMPUTER  SIMULATIONS  AND  ANALYSIS 

Using  the  techniques  described  in  Chapter  3,  the 
computer  simulation  for  this  configuration  has  been 
carried  out  for  a  device  with  the  parameters  given  below: 
doping  concentration  Nfl  =  1x10 17  cm-3 
built-in  potential  VQ  =  0.72  V 

diode  height  lx  =  0.1064  ym 

overlap  length  •  1Q  =  0.0266,  0.02,  0.0133  ym 

diode  area  11=  1x10 ~5  cm 

y  z 

.  .  2  -1  -1 
electron  mobility  y  =  5000  cm  V  sec 

n 

overlap  ratio  1Q/1  =  1//5 

depletion  length  WQ  =  0.0998  ym 

The  configuration  used  in  simulation  is  shown  in  Figure 
5.2.  Here,  N^  is  the  doping  concentration  for  the  bulk 
semiconductor,  and  1  is  the  thickness  of  the  bulk  region. 
The  two  Schottky  metals  overlap  by  1q.  In  order  to 
simplify  the  simulation,  the  Einstein  relationship  is 
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assumed.  Because  of  the  antisymmetric  geometry,  it  is 
necessary  to  analyze  only  one  polarity  of  applied  bias 
voltage.  In  the  z-direction  the  diode  structure  is 
assumed  uniform,  therefore  only  a  two-dimensional 
simulation  is  performed. 

Figures  5.3  and  5.4  show  the  carrier  concentration  and 
potential  contours,  respectively,  for  a  typical 
antisymmetric  diode.  With  zero  bias,  i.e.  graph  (b)  of 
Figures  5.3  and  5.4,  the  geometric  symmetry  of  the  diode 
is  reflected  in  the  contours.  The  electron  concentration 
is  normalized  to  the  doping  concentration  Nd.  Areas  where 
the  contour  values  are  less  than  n/Nd  =  0.001  are 
considered  to  be  depleted.  The  area  next  to  each  Schottky 
metal  contact  is  depleted.  This  phenomenon  is  similar  to 
that  observed  in  the  one-dimensional  case.  Assuming  that 
the  length  of  the  depletion  region  surrounding  each  metal 
contact  is  approximately  equal  to  one  depletion  length,  if 
the  thickness  of  the  bulk  region,  i.e.  the  distance 
between  the  two  Schottky  metal  contacts,  is  of  the  order 
of  two  depletion  lengths,  then  it  is  expected  that  there 
is  no  undepleted  channel  between  the  two  ohmic  contacts 
under  zero  bias  (See  Figure  5.3  (b)).  When  a  small  bias 
is  applied  (Figure  5.3  (a)),  since  the  left  hand  Schottky 
contact  is  forward  biased,  the  corresponding  depletion 
region  is  retracted  toward  the  metal.  At  the  same  time, 
the  other  Schottky  contact  is  reverse  biased,  and  so  the 


Figure  5.3  Typical  electron  concentration  contours  for  the 
antisymmetric  diode,  (a)  V  *  0.15  volt,  n 
increases  in  the  central  region  (b)  V  *  0,  the 
contours  are  symmetrical  (c)  V  *  0.30  volt, 
electron  conduction  through  the  central  region 
turns  the  diode 'on. 
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corresponding  depletion  region  extends  further  out  into 
the  bulk  material.  Comparing  this  with  the  zero  bias 
case,  it  is  observed  that  in  the  central  region,  the 
electron  concentrations  increase.  This  phenomenon  is  more 
clearly  observed  in  Figure  5.3  (c)  which  shows  an  exposed 
ohmic  channel  for  electron  conduction.  In  fact,  the  diode 
is  already  turned  on.  Notice  that  this  occurs  for  a  bias 
voltage  well  below  that  needed  to  turn  on  the  SBD. 
Instead  of  carrier  conduction  through  the  metal  contact, 
this  new  conductivity  modulation  device  provides  a  voltage 
variable  ohmic  channel  for  conduction.  It  will  be  shown 
later  that  the  I-V  relation  is  exponential  for  both 
polarities  of  applied  bias. 

Figure  5.4  shows  the  corresponding  potential  contours 
for  reference.  Similar  to  the  Gap  diode,  because  of  the 
uncertainty  of  boundary  values,  it  is  difficult  to 
accurately  describe  areas  where  the  Schottky  metal 
contacts  meet  the  ohmic  ones.  Away  from  these  areas, 
along  the  metal  contacts,  the  potential  becomes  more  and 
more  uniform.  From  these  contours,  the  electric  field  can 
be  found.  In  order  to  obtain  an  accurate  picture  of  the 
current  density  distribution  not  only  the  electron  contour 
plot  but  also  the  potential  contour  plot  muat  be  taaer 
into  account  (See  equation  3.2). 

Figure  5.5  shows  I-V  curves  for  different  over 
the  Schottky  metal.  Also  shown,  for  compar.s 
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I-V  curve  for  a  conventional  back-to-back  SBD  pair  with  V. 
=  0.72  V.  This  graph  shows  that  the  turn-on  voltage 

responds  sensitively  to  the  change  of  the  overlap.  More 
overlap  results  in  a  higher  tum-on  voltage.  Figure  5.6 
shows  the  corresponding  curves  on  a  log  scale.  Straight 
lines  in  this  graph  indicate  that  the  I-V  curve  of  the 
device  is  exponential,  independent  of  overlap.  Moreover, 
the  slopes  differ  from  the  back-to-back  SBD  pair.  With 
the  thickness  of  epi-layer  as  a  parameter,  the  resulting 
I-V  curves  are  shown  in  Figure  5.7.  For  a  thin  device, 
the  channel  is  tightened  up  and  electron  conduction  is 
restricted,  resulting  in  higher  tum-on  voltage.  The 
corresponding  I-V  curves  in  log  scale  are  shown  in  Figure 
5.8.  Again,  all  are  exponential  vut  with  different 
slopes . 

Figure  5.9  compares  calculated  capacitances  for  a 
typical  antisymmetric  diode  and  a  corresponding 
antiparallel  diode  pair.  The  capacitance  of  the 
antisymmetric  diode  is  much  lower,  since  the  current 
conduction  is  due  to  conductivity  modulation  rather  than 
to  the  Schottky  rectifying  effect. 

5.3  CONCLUSIONS  AND  COMMENTS 

According  to  the  simulation  it  can  be  concluded  that, 
for  the  antisymmetric  diode, 

1.  The  I-V  curve  is  antisymmetric  and  exponential. 

2.  The  tum-on  voltage  is  adjustable  by  two  parameters: 


miiiiiirjiirniirmiii'.iiiiimiiiiiiiiiiimiiiiiiiiiiii 

liiHiit^iriii^iiiiiriiiMiiiiniiiimiiimiiinniuil 

......  ..  .  i..~/ !..■■■■■ 

I  ■■■■■■ Jlk/ll  •' i  i ■■■■■■■■ !■■■■■■ ■■■■(!  !■■)■ ■■■■■■■■■■■■«■■ 

iaiiiiiiviii";r.iiiiiiiiriiiiii|iaiiiiiiiiiiiiiiiiaiimimia 

iiir.iirjii.jiiiiiiiijiiiiiiiiiiiiiiiiiiiimiiiiiiiiiimi 

I  ■  4  m  m  m  mmm  w  ■■■■■■  m  m  m  m  m  a  ■■■■■■  ■■■  ■■■■■■■■■■■■■■■■■■■  m  m  ataaia  ■■ 

i «■■■* ■■I*! ■■■■■■■■■? ■■■■!Saa 3S ■■■■■■■■■■■■■*■■■■■■■■■■■ 

iS'i!S'(iSSSSii!SSri!!S!SSSiS!iiiiSiiSSSiSi!SiSiS!SSSSSSSiu] 
lllll!illllllllllll/illllllllllllllllllllllllllllllllllllllll| 


liiiiiiiiiiiiiiiijiiiiiiiinm  iiiiiiiiiiiiiiiiiiiiiiiifiAM 

liiiiiiiliiiiiiir.iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiil 

llllllllllllllflllllllllllllllllllllllllllllllllllllllllllll 


iiiiiiiiiiiiniiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiim 

lllillllllKIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIilllllllllllllll 


iiimiiijiiiiiiiimiimiiiiiiiiiiiiiiimiiiiiiiiiiiiiiii 

lllllllllllllllllllllillllllllllllllllllllllllllllllllllllll 


imiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiimiiiiii 

llllllllllllllllllllllllllllllllllllllllllllllltllllllllllll 


**  9  r  ■###&&* _ ##&&&&&*’  y&60g&fr* 


1  i-Kyi”, 


-ku*  *  '  »-  »-  *-  *. 


Figure  5.7  Computer  results  of  I-V  curves  with  the  epi-leyer 
thickness,  1^,  as  the  control  parweter.  All  other 

paraaeters  stay  the  sane  as  before. 
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(a)  the  overlap  of  Schottky  metals. 

(b)  the  thickness  of  the  bulk  region. 

3.  Low  capacitance  can  be  expected .  other  properties, 
such  as  temperature  sensitivity  and  burn-out 
behavior,  should  be  similar  to  those  experienced 
with  the  Gap  diode. 

From  above,  the  new  diode  should  be  useful  as  a 
subharmonic  mixer  and/or  limiter. 

From  the  standpoint  of  design,  the  thickness  of  the 
bulk  semiconductor  should  be  no  more  than  3  depletion 
lengths.  The  reasoning  is  similar  to  the  gap  length 
restriction  for  Gap  diode.  Furthermore,  the  overlap  of 
Schottky  metals  should  not  be  taken  to  be  too  large,  since 
this  will  produce  a  low  shunt  resistance,  and  the 
resulting  device  will  be  a  resistor.  On  the  other  hand, 
too  narrow  a  channel  will  increase  the  turn-on  voltage 
and/or  capacitance,  and  the  advantage  of  this  new  device 
will  be  lost. 
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6.  CONCLUSIONS 


The  purpose  of  this  work  is  to  develop  a 
two-dimensional  simulation  program  capable  of  modeling  the 
static  electrical  characteristics  of  semiconductor  devices 
with  complex  geometries  and  made  of  arbitrary  material. 
The  simulation  is  intended  for  operation  on  small 
(desktop)  computers .  Another  purpose  is  to  use  the 
simulation  in  the  investigation  of  some  new  microwave 
device  designs. 

6.1  COMPUTER  SIMULATIONS 

The  program  that  has  been  developed  has  the  following 
properties: 

6.1.1  User-Oriented  program 

Information  is  entered  for  each  node  of  a  mesh  which 
describes  the  device  to  be  analysed  by  the  finite 
difference  method.  The  information  to  be  entered  includes 
the  neighborhood  node  identification  numbers,  the 
distances  between  nodes,  and  the  status  of  the  node  (i.e. 
is  it  on  a  boundary,  in  the  interior  of  the  bulk,  or  on  a 
symmetry  line,  is  it  Schottky  contact  or  ohmic  contact,  is 
it  forward  or  reverse  bias,  etc.).  The  carrier 
concentration,  the  potential,  and  the  current  density  at 
each  node  interior  to  the  device  can  be  found  as  well  as 
the  total  current  for  each  applied  bias  voltage. 


6.1.2  Modified  Decoupled  Method 

Because  of  the  limited  memory  available  in  a  small 
computer,  efficient  use  of  this  memory  is  a  primary 
concern.  A  modified  decoupled  solution  method  has  been 
developed  using  a  two  stage-iteration  scheme  which  not 
only  saves  three  quarters  of  the  memory  but  also  yields  a 
converged  solution  of  similar  quality  to  that  obtained 
from  the  traditional  coupled  method. 

6.1.3  Limitations 

The  maximum  number  of  nodes  available  for  simulation 
is  limited  by  the  computer  memory  size.  Also,  the  allowed 
maximum  mesh  spacing  is  dominated  by  the  Debye  length. 
Most  of  the  time  consumed  in  obtaining  a  solution  is  used 
in  solving  the  inverse  matrix.  Using  a  small  computer,  it 
might  take  many  hours  to  obtain  one  I-V  curve.  However, 
the  trade-off  is  the  low  cost  of  simulation.  Since,  for  a 
given  barrier  height,  the  ratio  between  depletion  length 
and  Debye  length  is  a  constant,  from  the  argument  above, 
the  total  area  which  can  be  simulated  is  limited  by  the 
computer  memory  size.  However,  usually  the  area  needed 
for  two-dimensional  modeling  is  limited;  other  parts  of 
the  device  can  be  treated  by  a  one-dimensional 
approximation.  Hence  it  can  be  concluded  that  by  suitable 
manipulation,  it  is  possible  for  the  designer  to  simulate 
a  variety  of  semiconductor  device  structures  with  this 
program. 


6.1.4  Finite  Difference  And  Finite  Element  Methods 

The  finite  difference  method  is  simple  and  reliable  in 
programming,  and  much  more  effort  can  be  applied  to  the 
device  simulation  itself  instead  of  to  debugging  the 
program.  Higher  order  finite  element  methods  might  be 
used  to  improve  the  results  after  a  basic  understanding  of 
the  device  under  consideration  has  been  made  using  the 
finite  difference  method. 

6.2  APPLICATIONS  IN  TWO-DIMENSIONAL  DEVICES 

By  changing  the  geometry  of  the  metal-semiconductor 
diode,  we  have  shown  that  it  is  possible  to  construct  a 
family  of  new  devices  which  can  be  used  in  microwave  and 
mm-wave  applications.  Bised  on  the  simulations,  some  of 
the  characteristics  c  i  two  new  devices  have  been 
predicted. 

6.3  GAP  CONTROLLED  METAL- SEMI CONDUCTOR  DIODE 

4 

In  this  new  design,  a  simple  and  effective  means  by 
which  to  control  the  turn-on  voltage  of  a  semiconductor 
diode  has  been  found.  This  new  contact  technique  permits 
the  use  of  conventional  Schottky  metals  and  any 
semiconductor  material.  This  design  may  also  improve  the 
temperature  behavior  and  burn-out  rating  of  these  diodes. 
Furthermore,  in  fabrication,  this  diode  only  takes  a  few 
more  steps  than  the  conventional  SBD,  so  that  the  added 
cost  may  not  be  significant.  An  actual  device  has  been 
demonstrated  in  GaAs  '37). 


6.4  ANTISYMMETRIC  METAL- SEMI  CONDUCTOR  DIODE 

Another  new  structure,  the  antisymmetric  diode,  was 
proposed  for  microwave  and/or  mm-wave  applications.  The 
potential  advantages  of  this  device  include  adjustable 
turn-on  voltage,  balance  in  both  polarities,  improved 
bum-out  rating,  low  capacitance  and  low  temperature 
sensitivity  as  a  mixer  or  limiter. 

6.5  PROBLEMS  TO  BE  SOLVED 

Two  new  types  of  microwave  and  mm-wave  diodes  have 
been  modeled  with  the  technique  described  in  the  study 
presented  here.  Other  configurations  also  can  be 
investigated.  For  example,  a  three  terminal  device  is 
possible  in  which  the  gap  technique  is  used  to  obtain  a 
vertical  FET.  A  sketch  of  one  such  design  is  shown  in 
Figure  6.1.  This  is  a  planar- type  vertical  FET. 
O  .paring  this  configuration  with  that  of  the  Gap  diode, 
as  shown  in  Figure  4.1,  the  only  difference  is  that  on  the 
top  portion  of  this  device  the  Schottky  metals  are 
separately  contacted  to  become  the  gates  (G).  The  drain 
(D)  is  formed  by  an  ohmic  contact.  Apparently,  the 
gate-to-gate  length  and  the  epi-layer  thickness  depend  on 
V. .  and  can  be  calculated  if  the  gap  technique  is 
applied.  Because  the  drain-source  distance  is  much  less 
than  for  a  conventional  FET,  good  performance  in  high 
speed  and/or  high  frequency  applications  can  be  expected 
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Figure  6.1  Cross-section  of  planar-type  vertical  FET. 

The  gate  (G),  drain  (D)  and  source  (S)  are  of 
Schottky  contact,  ohmic  contact  and  N+  substrate 
materials  respectively.  Gates  are  reverse  biased 
with  respect  to  the  drain  (Vqj<  0).  The  gate-to- 

gate  length  and  epi-layer  thickness  depend  on  VGD 
The  gap  technique  can  be  applied  to  decide  the 
device  parameters  (refer  to  Figure  4.1) 


In  order  to  improve  the  simulation  work,  other 
parameters  such  as  temperature  should  be  considered. 
Refinements  of  the  V(E)  and  0(E)  relations  used  also  would 
be  helpful  in  obtaining  more  reliable  results. 
Furthermore,  the  properties  of  the  contacts  should  be 
emphasized  since  this  is  an  important  factor  for  achieving 
accurate  predictions  of  device  behavior. 

With  the  simple  and  low  cost  technique  presented  here, 
the  simulation  of  semiconductor  device  static 
characteristics  is  accessible  to  laboratories  with  even 
modest  computing  capability.  Simulations  are  a  powerful 
tool  to  aid  in  the  understanding  of  complex  device 
phenomena.  The  knowledge  so  obtained  can  lead  to 
improvements  in  existing  devices  and  the  exploration  of 
new  device  concepts. 

Finally  in  this  study,  it  is  assumed  that  the  third 
dimension  of  the  device  .under  investigation  is  uniform. 
However,  for  some  small  device  in  which  the  third 
dimension  is  of  the  order  of  the  other  two  dimensions,  or 
for  which  more  accurate  results  are  required,  it  may 
become  necessary  to  develop  three  dimensional  analysis 
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Program  1  t  Nod#  Information 

•>  tntril  Pr«fria  far  natal-Seaieeadaetar  Blade 
la  frlva  •,797if llaa  Caaf«,CU».ej».Cai<i-i«ilaW 
ta  *atteal*elilapet*a 
3a  *ahee8*e8aealealattea.9taee  ataalaa  ta  aaay  aalpr agraaa . tat 
4a  'the  llaa#  aftea  dakmtaf.laMlMi|  aa  a-M'i 
Sa  *e-9C»e9C»a-9r»e9».e-98are98a.e-prtreprtaaaBprepraae»» 

•a  *Clapaatere  lapat  far  Best  aalaalatlaa* a  » 

7a  ‘CVNaatara  aatpat  far  aaat  ealealatlea'a 
•a  *CVW77aatara  aalpat  far  reeerd.Vt9I-V.N(91»J*a 
•a  *Clap??aatere  lapat  far  reeerd*a 
19a  *Cae da. atara  lapat  aa9e  lafaraatlaa*. 

11a  *Caede7aperaeeeet  atara  lapat  aada  lafaraatlaa* a 

19a  *Caada-a t— parary  atara  aada  lafaraatlaa* a 

13a  *CJ#aaa4ar  far  aaarraat  aalaalatlaa*  a 

14a  ‘Cfafar  a -AT  atara  taaaparary  reeel ta* a 

13.  'Hadee-laalda  a  tart  fraa  at  a  ad  laaraan'a 

13a  *  Nedee-fcaaadary  a  tart  fraa'.r-l  a  ad  4aaraaaa*a 

17a  *raad  a ad  prlat  a tart  fraa  laft  ta  rlfht.tep  ta  kattaa* a 

19a  dap  •aaraallaad  da  a  ad  dy*l'i«tt  3999 

13a  dap  *aa  aada  aa  aaa  alda.laaart  rare*  await  3999 

<9 a  dap  *la  aaaa  af  syaaatrle*i«ett  3999 

tii  dap  *alaaya  aala  it  9-8  ayaaaatrle'iaaalt  3999 

It a  dap  *aad  aaa  aerth  alda  far  aalaalatlaa* iaeit  3999 

X3a  fad  9|dla  94131 

t4a  aat  *tatal#  af  aadaa-*.T 

ZSa  aat  *dx<*.r,aat  ‘fya,llfila-J>P|«alr-})9 
tt.  fat  alt.ald.I.at 
97a  art  714,*da-*,P.*ee* 

93a  art  714,*dy*.9.*aa* 

23a  eat  'mm*,  aadaa  la  V  far  aalaalatlaa* ,Y 
39a  dta  X(Vl|dUa  BIV*tl 

31a  dla  Ztlllidla  V4tl ladle  J(8a41idle  9(41 ,dla  X9I31 

39a  -  .791«iS999>3aaat  •dapta9-*.H|1.9a-19/<19.3at>. 99a-14)aH>N 

33a  aat  'Cftaapa  file  Caada7*  .99 

34a  If  aap<99}-*VU* igte  *ahaa§e  aada* 

33a  aat  *S-ayaaatrla  er  R-aaaaayaaalr la*  .za 
38a  rraad  l.liaprt  l.Zt.'aed* 

37a  «ek  'title* 

33a  dap  'left  ta  rlfkt.tap  ta  kattaa*  t watt  9999 
33a  *aat*aaat  *aada»*,d«91 
49a  If  dl91<9|Jf91N.t*9*m 
41a  If  Jt#l>4|  Jt#J  Wta*I*>V* 

49a  fad  4«lf  dlt»9idap  'aada  laatda  r'MtrdOiaalt  3999apta 

43a  fad  lidap  *aada-9aaadary  r*»atrCL> iaait  3999 

44a  far  K-l  ta  4 

43.  BKItKlMHKl 

48a  aarat  K 

47a  yta  *9.8* 

48a  M.S'alf  M-i  er  H-taJap  3 

43a  aat  *tM9M  aadaa  aklft*17VK9  ar  98* .99 

39a  If  eap<94»-*VKS*ifte  •■91ft* 

31a  aat  *Caat  aada  «*.J(1I 

39a  aat  *Maat  aada  **.4(91 

33a  aat  *3aat9  aada  f*,JI3J 

34a  aat  *Nart9  aada  #*,JI41 

38a  *9*  aaat  *91 a tea aaa  9aaa  9a  taat  9aaT* ,94 

38a  If  aap<94>**VK8* i»ta  *3* 

37a  aat  *da<t>-* .Bill 
38a  aat  *daOO«*.D(91 
83a  aat  *dy(81a* .8131 
89a  aat  *dy(Nl-* ,BI4ia9te  *3* 

81a  *9.S*aaat  ■<8>eftettly?<BlheleV.I9 
89a  aat  *191 . <Bl-»iaaed7<9lr Beaded* .99191 
93a  *3*atf  aap(Z4)«*N*|X4I1.91»ltajap  3 
84a  aat  *aa  ayaaalrla  aat at  VKf  ar  MB* .99 
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«i  »f  M»(M><>vn>iIv<>ni3i 

99a  fat  M.I.Iiarl  714, J(f] 

97a  if  Wi.nort  714igts  >3.9* 

fat  Ii,4f4.l,li«rt  714,J(l).Jt3),.)l3),J(4J 
99a  art  714.0(11,9(31. 0(31.014) 

73  «  *3. 9*  •  fat  itSiurl  7I4,ufUI> 

71a  aat  ‘ItfaraatlM  Is  CsrraatTVCO.HO* ,A3 

7Z>  If  Mg(M>«aV(laift*  •arraf1 

73a  *Capy  Ns4s  lafaratlw*  a 

74a  aprt  t,YI,JI*l,Mll,M,anfa 

73a  *  *>X3|lf  IK*T||la  •prl* 

73a  |li  *aat* 

77a  *prt*idap  >1  all  aa4a  Infamatinn*  inait  3333 
73a  rraaal  1 ,1 laraad  l.Zliaa  and  I, *399* 

73a  gsh  •tills* 

33a  ataaimf  l,VS.J(a),0(a),xa 
31a  fat  f4.3, start  714, J(3) 

33a  if  JI3K3|fat  a 73 1 art  714,cap(X3>igts  *9* 

33a  fat  3s. 4f4.3. start  714.J(1K.J(3),J(3),J(4) 

34a  fat  4s.4f3.3.start  714,0(11,0(31,0(31,0(4) 

39a  fat  st4tart  714,aap<X3> 

33a  ft*  *3* 

37a  *399* aaat  Ms  ail  Isfsraatlsa  esrrast?*,A3 
33a  If  nag<A3J-*Vt»*,gtn  *3.3* 

33a  f la  “ahaii ga  aa4a* 

33a  *3.3* ifnd  Stdap  'Mark  lag  tha  landing  Nsdas*  I  as  It  3333 
Ota  far  1-1  ta  V 

33a  dap  '(.fading  Nadar  as  LlaarMstrOliaaait  3333 
33a  aat  3(11 

34a  dap  *Nada  Asa It  Isa*  aa  Llasr*3atr(|)|«alt  3333 
33a  sat  Kill 
33a  aaat  I 

37a  rraad  l.ltaraad  1,33 
33a  aa  sad  1,*3* 

93a  far  f-i  ta  V 
133a  far  J-l  ta  Sill 
131a  fat  a4,S|«art  714, •• 

133a  Mil  J 

133a  If  f Iglifta  *3* 

134a  *4*asraad  l,VS,JtB),D(B),X3 
133a  If  J(31-3(l*lltaf9  lifat  Mart  7t4(aast  I 
133a  *3* a fat  f4.3, start  714,J(31|Sfg  1 
137a  fta  *4* 

133a  aaat  I 

133a  *3* a  fat  /tart  714 

113a  aat  ‘Anything  Mraag  la  3rlataat7  VII .M3*. A3 
lilt  If  sap<AS>«*VtS*tgta  *3.3* 

113a  rraad  3,ltaprt  t,3.9,V,T,Xla).Atal,L,Ha9.tf.N.H.V,*aad* 
113a  fat  a93|«rt  7l4,*Nadas  far  Carraat  la  and  Oat* 

114a  fat  3,f3.3.4s,f3.3,s4,f3.4,15a,f9.3,4s,f3.3,a4,f3.4 
113a  rraad  3,1 
113a  far  1-1  ta  V 

117a  *agala*a4lap  *Llas«*Astr(|}|salt  3333 
113a  aat  Mat  aada  an  iapat*,A 
119a  ant  *Sad  aada  aa  lapat*,0  . 

133a  aat  »aar anil sad  dlataasa.rlaO  an  lapat*,C 

131a  amt  *lat  aada  an  aatpat*,0 

133 a  aat  *3ad  aada  aa  aatpat* ,C 

133a  aat  'aaraallsad  dlstaaaa,rl33  aa  aatpat* ,F 

134a  art  714.3.3,9.*  d«* ,3,0,1,*  d-*,f 

133a  aat  *3aaa thing  MraagTVIS.M*  ,A3 

133a  If  sap(A3>-*VC3*|fat  sS3|«rt  714, Mast  lias  la  nraag*igta 
137a  aprt  3,A.0,C,0,C,F.*aad* 

139a  aaat  t 

133a  aat  *9ld  yaa  ahaaga'Caadaf*?* ,A3 

133a  If  aap(A3>«*Vt9*tdap  *Chawga*Cnada7#* taalt  3333ihaaptatp 
131a  *a*aaapy  *Caada*,*aada33* 

133a  gat  *ahao3*.3,9 
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ml  ••ktft  Jlll,»lilf*,M 

If  aa*<04>r*VCi*i«ta 

far  1-t  ta  4 

JIIHMIII 

If  JIlMdMIli 

•aat  I 

•ta 

•Iff. far  I-t  ta  4 
aat  J(ll>a»t  0(11 
aast  I 

•«•»•  •>*• 
far  1-1  ta  3 
aat  Mill 

Mat  I 

Mr!  4,VI,JUIlMallM,W 
aat  *aaa  aara  Mat  ta  ItfMI 
If  aaa<A4>a*VCS*ista  *10* 

•ta  *10* 

•lr* aaaay  *Caa4a-*.*Caada* 
rat 

*aart*« 

rraaO  t.lirraa*  4,liaraa4  l.Ztiaart  4lZll‘aaf* 
aa  aaf  1 , “pr* 

* •#*••>*  "aarraat  aa4ar>?,tf  aa  aara  aat  •* , J 
*al*<araa4  I.YI.Jlal.MU.M 

If  J»J(01ia#rt  4,V*,J(a|,Diai.*«,*aa4*i*ta  >1* 
aat  *m4M  la  araaa?*,ft4 
If  aap( A4> •  •  VH* i aat  JI4I 
•at  'Jill  la  «ra^V,M 

If  •••€«•) --V(l*iaat  Jill, aat  Jltliaat  J(3Itaat  JI4) 

•at  *0(91  la  ■rn(f',M 

If  HflM>*‘Ktaiaal  Odliaat  OlCliaat  0(31»aat  0141 
•at  *>•  la  araag?* ,11 

If  aa»(M)>*V(S*l*  *>Miaat  Xtdliaat  Mltlimt  M(» 
•art  4.V«,J(a),ft(ai,M,*aa4* 

•ta  *M« 

*pr*iaaay  *Caa4a-*,*Caa4a* 
rat 


Program  2  *  Print  Nods  Information 

4a  file*  Cm* 

la  Ola  Z«til|4la  44131 

* a  Ota  V*III 

3a  Ota  44131 

4a  4ta  JCIa4l 

Sa  4la  1(41 

•a  Vl'i4a»  >1  alt  aa*  lalarMliw'iwU  MM 
7a  rra aO  l.lta r— 4  l,Z4iaa  aa4  t,*39S* 

•a  fat  'Ulla' 

•a  ■fair— 4  t,V*,3(al.4tal,M4 
14a  fat  f4.4, start  714,3(41 

11a  tf  JI4K4|f»t  o7*iart  ri4,aatU«)||ta  *V 
ICa  fat  4a,4f4.4.start  714,3(11.31*1,3(31.3(41 
13a  fat  4s,4fS.t,siart  714,04 tl, 41*1 ,4(31, 4(41 
14a  fat  «*4|«rt  7t4,afU*> 

ISa  |t<a  *4* 

14a  •tttla’afat  «44alf  Mf<n>S‘4‘l«t  714,*4yaa»trta  4»>astry‘i|M  * 
17a  aart  714,‘Wsssy— trls  inastry* 

14a  fat  «7,staart  714,*aa0a  #* 

14a  fat  all. mart  714,'aaOa  la  C.M.4.H* 

Na  fat  a*4.mart  714,’aaraal tsaO  fiaUtaa* 

41a  fat  «34«art  714,*C4)akattly4(4}l»la?1 

44a  fat  «44i«art  7t4.»Cri.<*l-4taaaO»C*>raaaOaO»» 

43a  art  714. *aa  sWaastrla  astaT* 

44a  rat 


Program  3  1  Main  Program 

It 

«•  inn  eM4*,ei«,M,u*leM,eM,eHynlevM«,erlej* 

«.  rraa*  l.tlinM  t.P.I.V 
3.  «■  KVIrfla  HYMIinI  11m. 

«•  «la  MI3l|leU»Mt*«|  « 

3t  m<4  t>TlH*l,IIH,LAM.IIi« 
t<  a#  *MJaat  twlM  V.n<  WiMlt  (HI 
7.  ml  Viaat  >«I,l|M.tcWUt.M.nml4)M 
3a  rraaC  t.laaprt  tl7ll,V,r,Xli)lMH,i,H,MiM>* 

•a  41a  VHIMm  jlta4lt«a  t|4l|«a  X3I3I 

III  Ma  VIlaHlaMa  NUaNlaUa  eiH)i«ia  FMliMa  ain.nl 
It.  If  VMirna  I,liaraa4  3,V|al,Ma| 

13a  't'lPraaa  l.liafal  ll7,l,V,T,Ill|ll|a),L,H,|,a,N,N,V,*ana< 

13a  rrma  7,tnprt  7,r,llV,Tll|iMli|,L,N,»,»,N,H,V,*ma* 

14a  1#  V*I,1M  V.Mal 
13a  rraa4  t,tiaraa4  l.Zt 
tta  m  n4  l,*!* 

17a  'I*>araa4  l,V4lJ(ai,||ala» 
tta  JIIIMilf  *>«||t#  •»• 
tta  xat i,t]>M(i/.atM>a 

tta  If  •an<A«)-*SP*iaaa<-V8>>nfAl|V-H>V(Alttf  V<4iam<V*»W*)|-«>VI»l 
tta  If  aaa<l4)«'tt' |4a|  ‘Hmg'ialp 

tta  If  a^CMFa^ff'ia^'MnmiliVmililf  V<l|l>M*l|»m*l 

33.  If  iaf(M)-aM‘|4aa  'Nraai'iata 

34.  If  aafMa)*'W* MI  Al  •  -31VI  A] » If  V<liaaa<W>M*l|-V-ll>V[l) 
tSa  If  aap<A«»>*33*|lMHAl|3>V(Alttf  V<«  laapCVClMI  A]  .-V1VI  A] 

tta  |la 

*7.  *t*alf  Vaitfta 

Ha  laa  A.C.P.Mata 

tti  rraaf  3,iiafrt  3,IHI,,a*4a 

33a  law  Ml 

31.  aal  M)F 

3Ca  far  1*1  ta  N 

33a  niimii 

34a  aaat  I 

33a  *M*  a  laa  I.Cirraal  3,l|araa4  3, Alai 

38a  far  l-t  ta  N 

37a  vtmeiu 

33a  aaat  I 

33a  aal  MIT 

43a  laa  A.Caaff  A.aAala  *#-AC*, 3.3 

41a  ara  F-Clf arraaC  3,laaprt  3.7lal.*aa3a 

43.  ahala  ■•-A#*. 3. 3 
43a  law  AM 

44.  rraa4  3.1»araaA  S.Ptal 
43a  *4Vt*aaat  ATM 

43a  3 1C 

47a  far  !•!  ta  H 
43.  e*C( I 1*313 
43a  aamt  I 

S3,  XlC/IOlC.fat  alS.alA.Saart  714, »#»!•* ,3 

31.  far  !•!  ta  N 

33a  Vlll-etlllVtll 

S3,  aaat  I 

34a  If  V-3i3ta  *3* 

33a  laa  A.AiaAata  *a-A3*,l3,8 

S3a  taw  AM»rraa4  4,lia»rt  4,Atai.*aa4* 

87a  aal  AC  IP 
33a  far  1*1  ta  N 
33a  ritlMK!! 

33a  aaat  I 

31.  rraaC  A.laacrt  3.Vtai,Ntal.»aa«* 

33  a  *3*  a  If  Ola-Safta  *33* 

33a  fat  «l7,fl3.3iart  714,*Applla«  Vallafw* .» 

34a  Mata  *a-prt*.t,3 


Program  4  t  Coefficient  Matrix  for  Poisson's  Equation 

••  i 

It  ppm4  1,1iipm4  1  >Z* 

ti  ••  m4  1,'tPit* 

3.  *alai*raa*  t,VMtil,M*l,» 

4.  If  JIIKIiiU  ••»• 

s<  all  'M'<jm,jiti,j»i,ji«i,mti.mti,Nui,«i4i,ju» 

e.  (It  *ai* 

?•  aarat*ilf  V#t  ar  fl|li«taii  ’chaal'.t.fl 
St  ahaia  'ikmlM,!* 

ti  aaaa.lf  plpX-ti-t/p7pSMtpt.p»)is>*  *a3a 

lit  If  p3p4«t|-l/p3p*M(pS.pt)ifak  *at* 

lit  If  plptp3p4rt|-t/p7pa-l/p3pSMCpS,pS) if  at  'H* 

tti  Sab 

13t  HCHIptlazpCV(p9)K)-l>»Ctp9)>Clp9) 
lft  rat 

is i  •as*tt/pe/«ps*pei>tii«  p2<fi-avcp2)>c(pti 

Itt  If  pl-ptt*MI  pt.pt)  MCpt.pt) 

17 t  If  pt>t  aat  pt*pt»ftM(pS,pt) 

Itt  */pt/(p3«pt>)t|tf  pKtiCCptl-SVtpllICCpt) 

It.  If  pl7titMIpt.pl  1 
tti  rat 

ti.  ‘arit/pe/cpityeniiif  p4<ticcpt)-«v(p4i>ccpti 
St.  if  pfMtMMpe.pti 

».  t/p7/Cp7*pt> >S|lf  p3<(|Clpt)-tVCp3)>C(pt) 
tf.  if  p3>tiMpe.psi*s>t(pe.p3i 
tSt  rat 
Mi  *a4attab 

17 1  |it 
Nt  rat 

Mt  aKa.if  r t»t 1 1/ . tSSSMC i rat 
Mi  SM7S 

ait  If  plp**ti <Vlpl  1-Vtpf  ))/(pS*pS7M 

3t t  If  p3p4rtt<VCp3)-V(p4»S<p7»pS>>B 

33t  \<M*SS»(itf  rt-titta  a\a 

S4.  CSe7S*l.la7«*t/3ttt“f»/Cl»f«/3»tt>“f>M 

3St  .  HSMtl/avlaMMH 

as.  WRWirat 

ait  'ViU  C<-attt|ta7/3tttMi.tt3M*t/3aacM>t 

as.  tf  (>aH|h?/(Ml.MSMt|/l»M«lat» 

at.  S/SMirat 


Prograa  5  t  Jacobian  Matrix 

i<  ■■-***. 

1.  If  \M||U 
t. 

3a  rraaf  I,l|Kwf  1,21 
4.  MMf  l.*#4* 

Sa  'M'lirMf  t,VI,J|ll,0|ll,ai 
It  If  JIIXIill* 

ft  «U  'P3*<Jlll.jm.J(31.J(4],PO(ll,PO(tl,18(31.10(41,Jlll7 

•a  |t«  •**• 

•a  •MSrraaf  S.tllfrl  S.tdl.'aaf* 

13a  If  Ml  flfttPpaa*  4,1  itpaaf  4,MH|Jip  3 

11a  laa  l,eiif|  3|«fg  tiatala  *a*Wa',M 

Ita  af|  3(lav  «M 

13a  rraaf  S.lirraa*  4.1 

14a  far  1*1  la  M 

19a  spaa*  S.Ftml 

18a  aat  tf>C 

17a  sppt  4,CIal,*aa** 

Ma  ant  I 

Ma  rraa*  4.1tsP*a*  4,M*I 
Ma  praaal  i.ltaraa*  1,2» 
tla  far  1-1  la  N 
tta  spaa*  l.vs. Jlal.Sla) ,Xl 
Ma  If  Jdlfltlsg  -I 

14a  all  *K*<Jlll,Jftl.JC31,Jt41,F8(tl,P8(tI.a8(31,iei41.Jtl]7 
Ma  far  J-t  la  H 

Ma  If  J-ltN<MJ.Il-MllK7aap(V(llK»MJ,ll 

*7 a  If  jaiiWMJ.llaaaCVfllK7MI3.il 

Ma  aaal  J 

Ma  aaal  I 

Ma  ppaa*  3.1 

31a  far  1-1  la  N 

Ma  spaa*  3.MH 

33a  far  J-l  la  It 

34a  MJ.llariJ17MJ.il 


Ma  aaal  J 
Ma  aaal  I 

37a  alaaia  “a* aat*. 1.43 

Ma  *#*ataa  tippaa*  l.lispaa*  1 ,21 

Ma  far  1-1  la  N 

41a  spaa*  1.V1, Jlal.llal.Kl 

41a  If  Jl 11*1  a Jap  -1 

41a  all  *K*<J«il,J<tl.Jt31.J(41,FMtf.MMI,M(3I.MI4l,Jlll7 
43a  Ml.H-MKaapCVIIlKI7Ml.il 
44a  aaal  I 

48a  aka la  ’sHaal* .1.43 
48a  *#3*afaP  l-t  la  4 
47a  I  *47  Jt  if  pi -Malta 
48a  if  plf-lapla  *••• 

48a  if  Pl-liCaspCCVfpl»/pjatMp81-Mpll>7Mp8,pll 
Ma  if  Pl-ligaa 
St  a  *M*aaaat  I 
Ma  rat 

S3a  *K*alf  Pl-4al/. IMSTKiPat 
Ma  17478 

Ma  If  plpt#l|(Vfpll-Vfpll>/(p8Pp8>74 
Ma  If  p3p4Pl|CVfp31-V(p411/<p7ap8>7a 
87a  1 <44*887 78|lf  pt-tafta  *\t* 

Ma  (8878*1.  l*7<K8S/3M8*477/<t*<8/38MI*47  74 
88a  .IMM't/34>IH4KM 
Ma  4/87K|pat 

81a  *\t*»lf  C<-3M4|M7/3M174| . 48884*8 /34r  14C4C78 
Ma  If  t>MMtta7/I74i .  —884 1/3*Plata7*178 
83a  4/878 iPal 
84a  *13»*a 


Program  6  :  Coefficient  Matrix  for  Current  Equation 

•  •  ‘e-MnM/.«aS/2>* 

li  rrtad  l,lnrH<t  1,Z» 

2.  aa  aad  1,'rat* 

3>  *al* .araad  1,V*. Jin), Dial, X* 

4.  If 

3.  ell  'u,Om,jm,J(»,J(41lrail)>»(2)l»(3],»(41,Jlll) 

S<  |la  •■»• 

7.  "rat*. if  flj3. chair  •i-MM,lt 
Si  chain  *ehaa2*,0,3S 

9.  •an*. 

10.  far  1-1  ta  4 
lit  if  pl«p9i9ta  *p* 

12.  if  pl'lijta  *p* 

13.  !*4>Ji9ab  *■* 

14.  if  pi<S|C(p9)*aH(pi>>eip9i 

13.  A(p9,p9J«9>A(p9,p9)  *. 

15.  if  X4-*»  ar  * i Jap  t 

17.  if  eap<Xf(3)>a*Va  aad  1*4  and  pl>4|-2a>a(p9,p4l igta  *p» 

IS.  if  pl>0»-9>ACp9.pll 
19.  *p*.aaat  I 
tt.  rat 

tl.  •■•.if  r0>tiaap(AVIpl))/pJ>9.rat 

tt.  abaCCVIpl !-V(p9))/pJ) IC.if  r0-2i3ta  *\* 

23.  CM7S.l.la7(Itt/3tta‘4»/(l*»/3(tU*4»F 

24.  . 0239F •2/3rrltri7>0 

23.  •0F*.0/pJaasp(F/0/ta<Vlp91«V(p|»>>9 
29.  rat 

27.  *\*.if  C<-3200|2a7/3200>F|.0239F*2/3nrlEFEF>D 
29.  if  E>3ttti2a7/C>F|.atSSF»t/3»rCl)ta7<‘2>9 
29.  sta  *0F* 


Program  7  t  Printout  Results 


•  ■  'fpri’iM  c3,el0.4,o2,4e,e3,el0.4,cf,el0.Z 

1.  «rt  N-‘,H 

2.  If  rf-f,fet  elf,f3.t»wrt  714,*aobtlity-*,f 

3i  f>t  c30i  If  714, ‘Einstein  tiutlm* 

4«  If  rl*li«rt  714,‘(«MrallMd  Nobility  and  VeUilty* 

St  if  rl*!i»rt  714,‘Thlo  Layer  Nobility  and  Velocity* 

St  fat  e28.elt.2,e3,if  rfrf,vrt  714, ‘Rales.  tlae-*,rl,*see* 

7.  toa  C,F)l>*>f iefy  1 
ft  rread  1,1, oread  l,Zf 
9,  *wl* toreed  1 , Yf , Jf SI ,D(o] , X* 

1ft  Jff)>C(Al,if  A<N,A»l>A,gte  M* 
lit  on  end  l,‘w3‘ 

12t  *«2‘, oread  l,ve,j(ai,oia],xt 
13a  J(91>FtBl,8-l>9,gto  ,»2‘ 

14 t  'oS'ifot  elfiert  714, *Vi-* 

ISt  1>« 

18.  rread  l,l,oreed  l,2f 
17.  for  1-1  to  V 

lft  90b  *opace* 

19.  90b  *VN* 

2ft  nest  I 

21 1  17A,sfg  l,fat  i net  714, •• 

22.  If  Zf— S*,fat  cttiort  714. ‘Last  line  to  Synaetrie  fsls* 

23.  fat  elf,«rt  714, ‘Hi-*, f>X 

24.  rread  1,1, oread  l.Zf 
23.  for  1-1  to  V 

26.  gob  ‘spoee* 

27.  9ob  *VN‘ 

28t  neat  I 
29 1  efy  1 

3f .  fat  ,nr  t  714,“ 

31t  If  Zf-*S‘,fat  off ,wrt  714. ‘Loot  Lina  ta  Syaaetrlc  Asia* 

32.  ‘Capae 1 tanee*  1 1 . 6e-19aM(T-X> IX 
33 t  fat  e2f ,el3.f,art  714, ‘total  ehnrye** ,X 
34 1  *Cerrant*,fat  /,wrt  714 
33 ■  rread  lf,l 

36 i  fat  elf,34s,e2f ,ert  714, "Cor rent  dans.  ln‘,*Cnrrent  dans,  ant* 

37.  fat  l,ell.3.42a,ell.3 

38.  tISIM 

39.  for  1-1  to  V 

49.  oread  If ,A,a,C,0,E.F 

41.  If  rf-fifte  *J1* 

42.  ell  MK*<A.8.FC.J1 

43.  ell  MK*<0,E.Pr,K> 

44.  gte  *SM‘ 

49.  *Jl*,esp((Vf  A1-VI91)/.  9299/21 /CFa<Mt8)-M(A))Hfa. t299at .fa- 191 J 
48.  eap«VtP>*V(E)}/.t2S9/21/FFsCNIEl-Nt0)}Hfa.tZ98al.8e-t9>K 

47.  *SH*.S«3>8,H*K>N 

48.  ert  714.1.J.K 
48.  nest  I 

Sf.  fat  e4f,art  714. ‘Average  Corrent  fans.* 

91,  If  eap<Ztl-‘N‘igte  *pl* 

92,  S-J/218,H-K/2>M,vrt  714.1,9/07-. S> ,H/(V-.S> 

S3,  <8*10/207-. 91  >Ntf],V>VIfl 

94,  gte  *pf* 

99,  *pl*,eet  714. 1 ,9/V.M/V 
98.  <S*M>/2V>NI91|V>Vttl 

97.  *p2*. chain  ‘eben2‘,f,89 

98,  ‘spaae*. 

98,  far  J-l  ta  Xtl) 

89,  fat  c8,t,art  714,“ 

81,  nast  J 

82,  rot 

83,  ‘VH*.tf  A>N,FIA-Nl>8,jap  2 

84,  CIA118 


xa82*s*3*asass23s3sja3:as2ss85  8a 


if  ■•«(l*lllfat  i«rt  714ir*t 
fat  f4.S.<tlf  Mi  f Ifliart  7I«,VII)||U  *T* 
if  rl«l  *r  l<l|t(ll«p(Vlll/.|Bt))Ji<irl  714,J|ft*  *T* 

Safe  '»• 

H(l]aap(KV(l)>> 7I4.J 
•T*.J»X>X,if  t'Tirat 
*VH* 

•JK*>afe»(tVlpl J-VIp*l>/p3>>r#|if  r4-2t»ta  *\1* 

<4279*1 .  U7(rf3/3IH*«»/ <1  •  (rt/IIM)*4)M 
.  (2SSr3*t/9al(rtr3)*»r4 

ap4ati.(ri*fe«tti7(rl/p4/tMVttll*VtfM»/f)i<m7t!-Mpi»>f4 
rat 

•\l*iif  r«*»H|N7/3ni>r3l.MStr3«t/»rl(rtr3)*nr4 
i f  r»3204 i2a7/r»r3i.  »»4r3*f/3*r 1*2*7*21*4 
ft*  *p4* 

•K'i 

•ir’itraaf  l.Yf , JC*I ,01*1 ,X* 
if  Jtll«li«U  *lr* 

ell  >M<Utll,Jltl,Jl»,Jt«l,mtl,PMt)ltMJI,«tl4l>JII» 
•2r*«r*t 
*Ni*>41Z>M 

if  plp2r4| (VI pi 1-Vt  p21 ) / <pS*p4)  Ht 
if  p3p4«t|(Vlpll*Vlp4l)/(t7*m>Z 
\  <HM*ZZ1 1C  t  if  r4*2ift* 
(W75»i.la7(ni/3tHM))/(l«(l/3tN>M»F 

F/OlKiret 

•\2*.if  C< *3244 1 2*7/3244 If  i  •42S9F*2/3*rlCFCF14 
if  03244)2*7/0 F • . 4299f*2/3fer laC* 7*270 
•  F/0>Kir*t 
•12947 
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